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In the name of the Lord of both wisdom and mind, 

To nothing sublimer can thought be applied, 

The Lord of whatever is named or assigned 

A place, the Sustainer of all and the Guide, 

The Lord of Saturn and the turning sky, 

Who causeth Venus, Sun, and Moon to shine, 

Who is above conception, name, or sign, 

The Artist of the heaven's jewellery! 

Him thou canst see not though thy sight thou strain, 

For thought itself will struggle to attain 

To One above all name and place in vain, 

Since mind and wisdom fail to penetrate 

Beyond our elements, but operate 

On matters that the senses render plain. 

None then can praise God as He is. Observe 

Thy duty: 'tis to gird thyself to serve. 

He weigheth mind and wisdom; should He be 

Encompassed by a thought that He hath weighed? 

Can He be praised by such machinery 

As this, with mind or soul or reason's aid? 

Confess His being but affirm no more, 

Adore Him and all other ways ignore, 

Observing His commands. Thy source of might 

Is knowledge: thus old hearts grow young again, 

But things above the Veil surpass in height 

All words: God's essence is beyond our ken. 



From Shahname 



By Hakim Abu'l-Qasim Ferdowsi Tusi (940-1020 CE) 
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Chapter 1 
Introduction 



It is not exactly clear when the very first idea of a computer was born: It 
seems that the first idea of an automatic computing machine was proposed by 
Charls Babbage in 1822. In the middle of 1930, Alan Turing by introducing 
the universal Turing machine, made it clear what characteristics computing 
machines needs to have, and how important the programming is. Shortly 
after his paper, the first computer was built. 

From 1947, onward transistors made the computer hardware develop very 
rapidly. So steadily the miniaturization pace that, in the middle of the sixties, 
Gordon Moore claimed that every two years the number of transistors used 
in a specific volume would doubl. This obviously had increased the speed of 
processing and the memory of the computers. This technology development 
went on until now, that we are not using the 30 tons computers anymore, but 
much smaller and lighter. However making computers and their hardware 
smaller, requires new considerations: when facing atomic scales, quantum 
effects play the major role instead of classical ones. Therefore if we want 
to have smaller computers, we might need to replace the classical computer 
technology by a quantum one. 

In 1982, Richard Feynman, mentioned that simulating a quantum system 
by a classical computer is very difficult Instead he suggested to use a 
computer which is ruled by Quantum Mechanics. Three years later, in 1985, 
David Deutsch proposed a model of quantum Turing machine [22] . This was 
one of the first attempts to reveal the extreme power of quantum computing. 
The interesting point is that quantum technology not only permits the use 
of very small size microchips, but it also gives the possibility to tackle cal- 
culations by a new generation of algorithms based on Quantum Mechanics. 
It was Peter Shor, in 1994, who made one of the biggest steps in this field 
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and showed how powerful such a quantum algorithm can be in factorizing a 
given number into its prime factors [52] . 

In order to understand the main idea behind such an algorithm, let us 
mention some simple facts. From a physical point of view, a bit is a system 
which can be in two different states; these states can be yes or no, right 
or wrong, or 1. A classical bit can be a mechanical on/off switch, or an 
electronic device capable of distinguishing a voltage difference. But if we 
choose an atom as a bit, in this case called qubit, then quantum mechanics 
tells us that the state of the atom can be a superposition the only two states 
available to a classical bit. Mathematically we write such a state as: 

\1>) = a\0)+/3\l), \a\ 2 + \(3\ 2 = l. 

Now, consider a classical memory of three bits: at a given time, it can be in 
only one of the eight possible states 000, 001, 010,. ..,111. On the other hand, 
a quantum memory of three qubits, can instead be in all the eight states at 
the same time:\ip) = a!|000)+a2|001) + - ■ - + 081111). Therefore if we increase 
the number of qubits, the information stored in them increases exponentially: 
a memory of n qubits stores 2 n classical binary strings simultaneously, while 
its classical partner does store only one string. This gives the possibility of 
performing the computation on the superposition of all the possible states in 
one step. This is called parallelism: a quantum computer can do a certain 
computation on 2 n numbers using n qubits in one operation, while a classical 
computer has to do the same computation 2 n times with different entries, 
or it has to have 2 n processors operating at the same time. This means 
that quantum computers can in line of principle perform exponentially faster 
than the classical ones. Of course, after a quantum computation of sort, 
the solution of a given problem will be contained in just one of the classical 
strings appearing in the quantum superposition; indeed, the very core of 
quantum algorithmics is not only to write a sequence of quantum gates in 
order to solve the problem, but also to devise protocols able to read the 
superoposition so that it gets projected onto the solution string with very 
high probability. 

However this is not the only spectacular feature of the quantum computers. 
The other aspect which makes them very much different from their classical 
counterparts is the quantum effect called entanglement. In 1935, Einstein, 
Podolski and Rosen, tried to explain the odd characters of the entangled 
quantum states, using the hidden variables [25] . It was however, John Bell in 
the sixties, the one who deeply investigated, the borders between Classical 
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and Quantum Mechanics. His celebrated inequalities [TO], that are to be 
satisfied by any local classical probabilistic theory, are violated by entangled 
states. That means that entanglement has no classically local counterpart. 

Consider the following state: 

M = ^(|01> - |10», (1.1) 

where Alice has the first qubit and Bob owns the second one. Here the states 
|0) and |1) are the eigenstates of the spin along the Z axis, with eigenvalues 
+ 1 and —1 respectively. If Alice does no measurements on her qubit, she can 
only guess that Bob's qubit is in the state |0) with probability | or, it is in 
the state |1) with the same probability. However if she does a measurement 
along the Z axis on her qubit, according to what she gets, she can say exactly 
in which state Bob's qubit is. For instance, if after measurement, she gets 
+1, then she can say that Bob's qubit is in state |1). While if she would have 
got —1, then Bob's state would have been |0). 

Now consider the state 

|^) = -J=(|0) + |1))®|0). (1.2) 

In this case Alice can say that Bob's qubit is in the state |0) without doing 
any measurement on her qubit. Indeed Bob's qubit state, before and after 
the measurement on Alice's qubit does not change. While in the first case 
the state of Bob's qubit was modified by Alice's measurement, in the second 
case it did not. 



This gives an idea that entangled states such as (1.1) are correlated in a 



way that the separable states as in (1.2) are not and that these correlations 
have no classical counterparts. 

As we mentioned before entanglement plays a very important rule in quan- 
tum computation and quantum information. Quantum dense coding and 
quantum cryptography are examples of that. However the most astonish- 
ing application of entanglement is quantum teleportation, which concerns 
sending a qubit from one party to another one. 

Suppose Alice wants to send the following state to Bob: 



\<f>)c = a\Q)+/3\l). 
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All Alice knows about the coefficients a and f3 is that \a\ 2 + \/3\ 2 = 1; that 
is the state \(f>)c is unknown to her. In order to send this state, they share 
one of the following maximally entangled states (later we will see that these 
states are called Bell states, and form and orthonormal basis in C 2 <S> C 2 ): 

I**) = ^(|00>±|11», 
I**) = ^(|01)±|10». 

Let Alice and Bob share the state |$ + ). Then the overall state would be: 
\4>)c ® \$)ab = ^(«|000) + a|011) + 01100) + /3|111» 

This can be re-written as: 

\4>)c®\®)ab = l(\® + }Ac®(a\0) + P\l))B + \$-)Ac®(a\0)-P\l))B 

+ \^ + )ac ® Hi) + P\0))b + \V-)ac ® (a|l) - P\0))b). 

Now if Alice does a measurement of the Bell basis on her first two qubits, 
according to whether she finds them in l^), l^) the state would collapse 
into one of the following states: 

|$ + )Ac®(a|0)+/3|l))B 
|$->AC®(a|0)-/9|l))fl 

\^ + ) AC (g)(a\l)+P\0))B 
\V-) AC ®(a\l)-P\0)) B . 

Therefore if she calls Bob through a classical channel (phone, e-mail, 
skype, etc.) and tell him what was the outcome of her measurement on Bell 
states basis, then Bob can do one of the following operations on his qubit in 
order to get it into the initial \4>)c state: 

• If Alice gets |$ + ), Bob does nothing. 

• If Alice gets \&~), Bob acts with <t 3 on his qubit. 

• If Alice gets |^/ + ), Bob acts with o\ on his qubit. 

• If Alice gets \^~), Bob acts with a^ai on his qubit. 

Where <j\ is the first Pauli matrix, and <t 3 is the third one. This way Alice has 
sent her unknown qubit to Bob. This process is called quantum teleportation: 
it amounts to sending one qubit of information using two classical bits and 
an entangled state as transmission channel. 
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Ever since, quantum entanglement has been one of the the main ingredients 
of quantum information theory, and the focus of many studies [37]. However 
very fundamental questions have still partial answers: when is it that a 
generic quantum state is entangled? How much entangled is it? Or, is its 
entanglement useful for quantum information processing or not? Indeed, just 
detecting entanglement in physical systems is one of the most challenging 
issues. 

Interestingly, beside rejuvenating the whole of Quantum Mechanics, quan- 
tum information has also stimulated a more abstract interest in an old math- 
ematical problem: the characterization of positive linear maps on (the states 
on) C* algebras of operators |36j . 




From a physical point of view, all actually occurring processes must be 
described by completely positive (CP) maps whose structure was fully char- 
acterized by Stinespring in the 50ties and whose central role in Quantum 
Mechanics was clarified by Kraus in the 70ies [4] . Their importance is due to 
the fact that they are much more than positive: let the system on which they 
act (A) be coupled to an ancillary system (B) completely inert and extend 
the action of the map from A to A+B by composing it with the identity ac- 
tion on B. Then, CP maps preserve positivity when they act on the states of 
the bipartite system A+B. Instead, only positive maps do not because there 
surely exist entangled states of the compound system A+B whose spectrum 
does not remain positive under a positive, but not completely positive map 
composed with the identity. 

Since one cannot physically exclude such statistical couplings, complete 
positivity is the only guarantee against the appearance of unphysical nega- 
tive probabilities in the spectrum of a transformed entangled state in the case 
of couplings to ancillas. The most renown amongst positive, but not com- 
pletely positive maps is matrix-transposition: as transposition, it preserves 
positivity, but fails to do so as partial-transposition, that is when acting on 
one party only of a bipartite compound system [36] . 




Thus, positive maps are not fully consistent as descriptions of actual phys- 
ical processes; however, exactly because they do not preserve the positivity 
of all entangled states, they can then be used as entanglement witnesses. 
Transposition is such a witness, which is exhaustive only in lower dimension; 
that is for two qubits (2-dim. system), or one qubit and one qutrit (three- 
dim, system). Otherwise, for instance already for two pairs A and B of two 
qubits each, there are entangled states of A+B which remain positive under 
partial-transposition: they are called PPT entangled states, namely positive 
under partial transposition, but entangled nevertheless. 
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The physical interest of PPT entangled states is that their entanglement 
content cannot be augmented by any distillation protocol operating on great 
numbers of them: they are also called bound entangled states as their entan- 
glement is somewhat locked in and cannot be extracted. The mathematical 
interest of these states is instead related to the notion of decomposable pos- 
itive maps. It was introduced by St0rmer and Choi in order to isolate a 
particular sub-class of positive maps and used by Woronowicz to character- 
ize all positive maps in low dimension (this makes transposition an exhaustive 
witness in such cases) [36J. 

Only a sub-class of higher dimensional positive maps is decomposable, 
that is can be split into the sum of a CP map and a CP map composed 
with transposition and these cannot detect the entanglement of PPT states. 
Therefore, one is interested in constructing new families of non-decomposable 
positive maps able to witness certain bound-entangled states. This activity 
has been pursued in recent years, among others, by the Horodeccy [36J, by 
Kossakowski and Chrusciriski (2TJJ, [21] in Poland. 

In this thesis work we have studied of 16 x 16 density matrices of two 
pairs of qubits, equidistributed over particularly symmetric orthogonal one- 
dimensional projections P a p,a,/3 = 0,1,2,4, for which Benatti et al. have 
been able to classify all PPT states, but only a few entangled ones among 
them. As PPT states are related to specific discrete geometric structures 
associated to a 16 point lattice, the issue is to understand whether bound 
entanglement might be related to a further specification of such structures 
and, if so, to characterize all PPT entangled lattice states. 

The organization of the thesis will be as follows: 

Chapter 2: We present an essential introduction of Functional Analysis 
as the basic mathematical tools dealing with Quantum Mechanics: Hilbert 
spaces and operators acting on them with a couple of most important norms 
on them will be defined. In this work, we shall consider only finite dimen- 
sional Hilbert spaces, H = C d ; therefore the algebra of bounded linear oper- 
ators defined on them will be isomorphic to the full matrix algebra Md(C). 

Chapter 3: We explain how quantum systems can be described using 
the mathematical tools of Chapter 2. Density matrices will be introduced 
as quantum states, which are, in general, mixtures of pure states with corre- 
sponding weights. Pure states, that is projectors onto Hilbert space vectors, 
are a particular case when the physical state of a quantum system is com- 
pletely known. 
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The von Neumann entropy as a counterpart of Shannon entropy will also 
be introduced; in compound systems, this quantity is used to compare the 
mixedness of density matrices of subsystems, called reduced density matrices, 
to the one of whole system. Throughout this work, we shall be interested in 
bipartite systems only. The definition of entangled and separable states in 
such systems will follow. We shall see, that as far as one is dealing with pure 
states, thanks to the Schmidt decomposition, distinguishing entangled and 
separable states is an easy task, which is not the case for mixed states. We 
will present Bell states and Werner states as examples of entangled states. 

Chapter 4: We explain what are the requirements for a linear map which 
is to describe a physical transformation. We start by giving some mathemat- 
ical definitions of positive maps and completely positive maps. Using the 
Kraus representation we will study the structure of CP maps. Then Choi 
matrix will be introduced, as it provides a technique based on its block pos- 
itivity or positivity to decide whether the corresponding map is positive or 
completely positive. We shall see that, though positive maps are not suitable 
candidates for describing physical transformations, nevertheless they can be 
used as mathematical tools to detect the entanglement. 

Many examples will be given throughout the chapter; particularly the 
Transposition Map, which will be used to define decomposable maps. We 
shall also see, how a generic positive map can be written in terms of the 
Trace Map and a suitable completely positive map. 

Chapter 5: This chapter focuses on variety of entanglement detecting 
methods. Partial Transposition, which is exhaustive in low dimensions, and 
PPTness as a necessary condition for separability will be presented. Using 
the Hahn-Banach separation theorem, we introduce entanglement witnesses 
and how they can be connected to positive maps through the Jamiolkowski 
isomorphism. The Reduction Criterion, which is based on the Reduction 
Map, already introduced in chapter 4 as a decomposable map, will be pre- 
sented. Further, we shall rapidly overview the Range Criterion, which was 
used to detect the first PPT entangled state in C 2 <g> C 4 , it will be followed 
by the notion of Unextendible Product Bases (UPB). We end this chapter 
by explaining the Realignment Method, which, unlike the other methods 
which manipulate only one subsystem, involves both subsystems to detect 
the entanglement. Many examples will be given throughout this chapter. 

Chapter 6: The last chapter contains new results. They concern a class 
of states over two parties consisting of n qubits each. The states studied are 
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diagonal in the basis generated by the action of tensor products of the form 
t 2 n<g>(jfl, ap = ®" =1 cr w , on the totally symmetric state ) G C 2 ™<g>C 2n . We 
first characterize the structure of positive maps detecting the entangled ones 
among them; we will show that their possible entanglement will be witnessed 
using a subclass of positive maps, namely diagonal positive maps. Then, 
the result will be illustrated by examining some entanglement witnesses for 
the case n — 2. Further, we will show how, for pairs of two qubits, being 
separable, entangled and bound entangled are state properties related to the 
geometric patterns of subsets of a 16 point square lattice. 



Chapter 2 
Preliminaries 



2.1 Banach and Hilbert spaces 

In this chapter we focus on the essential background of functional analysis [571 
[691 HI], an d review the basic mathematical tools which are used in Quantum 
Mechanics as we are interested to describe it using the vector and operators 
over vector spaces, in particular Hilbert spaces. 

Definition 1 Banach Space 

A Banach space is a vector space with a norm defined on it, with respect 
to which, it is complete. 

Definition 2 Hilbert Space 

A Hilbert space H is a Banach space with respect to a norm given by the 
scalar product (.|.). For every vector ijj G H, || ip ||= denotes its 

norm. 

We shall consider Hilbert spaces with countable orthonormal bases C 
H, the corresponding projectors Pi = satisfying the completeness re- 

lation: 52. Pi = i. 

Definition 3 Operator Norm 

The norm of an operator O on a Hilbert space H is defined as: 

|| O ||:= sup || 0\ip) || . (2.1) 
IMI=i 

It satisfies: 

|| Ot ||=|| o ||, || ] || = || O f . (2.2) 

O is said to be bounded if \\ O ||< oo. We denote the set of all bounded linear 
operators acting on H by B(M). 
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The linear combination and the product of two bounded operator is also 
bounded. Moreover all Cauchy type sequences with respect to the norm(|3|, 
converge to an element in B(M) [5] which makes this space complete. That 
means that B(M) is a C* — algebra. 

Definition 4 Spectrum 

The spectrum of a bounded operator O denoted by cr(0) ; is the set of all 
complex numbers X such that Aid — O has no inverse in B(M), where id is 
the identity operator. 

The projectors onto vectors of Hilbert spaces are the simplest examples 
of bounded operators on H with the spectrum a(Pj) = {0, 1}. 

In Quantum Mechanics positive linear operators are highly important, 
they are bounded operators with non-negative spectrum, indeed they can be 
considered as observables whose measured outcomes are positive. 

Definition 5 Positive Operator 

An operator O G B(M) acting on H ; is positive if (ip\Oip) > for all 

|V>)eH. 

B + (M) denotes the set of all bounded linear positive operators acting on 
a Hilbert space H: £+(H) C £(H). 

A set is said to be a cone if it is invariant to non-negative scaling. For 
any positive operator O G B + (M) it holds that aO G £> + (H), where a > 0, 
therefore the set B + (M) is a cone. 



2.2 Operators on a Finite Dimensional Hilbert 
Space 

Every finite Hilbert space H of dimension d is isomorphic to vector space 
of d-dimensional complex vectors C d . Fixing the orthonormal basis set {\i)}, 

( \ 
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i-th, 



any vector 



V j 

G H can be represented over the field of complex scalars: 

: G C d , 
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and linear operators are represented by d x d matrices. Then the physical 
observables are Hermitian d x d matrices with complex entries, and their 
algebra is denoted by M d (C), the entries of A G M d (C) being Ay = (i\A\j). 

Theorem 1 Spectral decomposition 

Any normal operator A, A^A = AA* , on C d can be diagonalized with 
respect to the orthonormal basis of its eigenvectors: 

A = J ^2a i \a i )(a i \, A\ai) = a^i) . 

i 

Therefore, A G M d (C) is positive semidefmite if all its eigenvalues are 
non-negative. Notice that matrix algebra M d (C) can be treated as a Hilbert 
space equipped with the Hilbert-Schmidt scalar product of two operators A 
and B in M d (C) which is defined to be : 

(A\B) = Ty{A^B), (2.3) 

where 

Ty{A^B) = ^2(i\A^B\i) 

i 

is the trace of A^B which is defined with respect to any orthonormal basis 
set {\i)}. The corresponding norm, called Hilbert-Schmidt norm: 



X>?, (2-4) 

can be expressed in terms of the eigenvalues af of the positive operator | A\ 2 = 
A^A. Another equivalent norm is the so called trace norm: 

d 

|| A \\i = Tr|A| = ^2 \ a i\- ( 2 -5) 
i=i 

Example 1 Trace map J&jj 

Let H = M d (C), where {F^f^ is a set of d x d matrices which are 
orthogonal with respect to the Hilbert-Schmidt scalar product, Ti{F^F k ) = 5^. 
They form an orthonormal basis set, so that: 

X = Y, Tr(XF})Fi, VX G M d {C). (2.6) 
i=i 



A || 2 = y/Tr(AtA) = ^ 
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Now consider the Trace Map, Tr d : M d (C) — > M d (C) defined to be: 

M d (C) 9lM Tr d [X] := Tr(X)id d . (2.7) 

Given an orthonormal basis {|a)}a = i, one constructs the so called unit ma- 
trices E a/ 3 := \cx)(/3\ in M d (C), which form an orthonormal basis in M d (C) 
with respect to the Hilbert- Schmidt scalar product. Therefore, using this or- 
thonormal set, the Trace Map reads: 



J2(a\X\a)J2\P)( 




a=l /3=1 




d 


d 


\fi)(a\X\ a ){fi\ 




a,f3=l 


a,(3=l 



since the two orthonormal basis {Fi}^ and {E a p}^ a /3=1 ^ can be transformed 

one into the other using a unitary matrix U E M d 2(C): E a/3 = X)j=i U a (3,jFj- 
The Trace Map can thus be re-written as: 

d 

Tr d [X] = ^2 E a/3 XE a/3 

a,/3=l 
d 2 d 

= E E U j,a3 U "/S,iFjXFi 

i j=l a,/3=l 
7,5 = 1 

d? 

j'=i 

Example 2 The simplest quantum systems are the 2-level systems, which 
are described by 2x2 matrices: A e M 2 (C). Most distinguished among them 
are the self-adjoint Pauli matrices: 

*-(;;). ^-(i-i). <-) 

wjt/i respect to the orthonormal basis |0), |1) o/ eigenvectors of a 3 : <7 3 |0) = 
|0), <x 3 |l) = |1), where |0) = (J) and |1) = (°). For £/ie ot/ier matrices we 
have <7i|0) = |1), <n|l) = |0), <r 2 |0) = and <t 2 |1) = i|0). 

The algebra of Pauli matrices is summerized by: 

CjO-k = o~jk^2 + iSjkiO'i (2-10) 
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where Sjki is the antisymmetric 3-tensor, 1 2 is the 2x2 identity matrix 
which we also denote by o~q. These four matrices are orthogonal in the sense 
of Hilbert- Schmidt scalar product, that is Tr(<Tj<Tj) = 5ij if i ^ j, when nor- 
malized, &i := they form an orthonormal basis set in M 2 (C). Therefore 
any matrix X e M 2 (C) can be written as: 

3 

X = ^(Tr(^X))^. 

i=0 

Example 3 Trace map in M 2 (C) 

Using the Pauli matrices as an orthonormal basis set we can write the 
Trace Map defined in (2.7) as: 

3 

Tr[X] = J2^Xa r (2.11) 

3=0 

In this case Pauli matrices permit us to simplify this map even more. Let us 
define the map S a : M 2 (C) — > M 2 (C) such that: 

S a [.} = a a [.]a a . (2.12) 

The Trace Map is re-written as: 

3 

Ti 2 [X] = J2Sa[X}. (2.13) 

a=0 



We shall refer to this map in the next chapters. 
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Chapter 3 

States and Operators 



In this chapter, we present some mathematical aspects concerning quan- 
tum states of multipartite systems. The presence of many parties introduces 
the notion of quantum correlations and entanglement. For the rest of this 
work, we consider only finite complex Hilbert spaces H = C d , where d denotes 
the dimension of the Hilbert space. States both in Classical and Quantum 
Mechanics are generic linear, positive, normalized functionals on the algebras 
of observables that fix their mean values. In Classical Mechanics one deals 
with probability distributions on phase-space; in Quantum Mechanics distri- 
butions are replaced by density matrices which extend the notion of quantum 
state beyond the single vectors in Hilbert spaces. 



3.1 Density Matrices 

Usually one is accustomed to states of quantum systems only as vectors in 
a suitable Hilbert space. However, more general ones need to be introduced. 
Consider an ensemble of projectors Pj := \ipi){ipi\ G -B(H) corresponding to 
the vector states G H, with weights pi. If the only knowledge about the 
physical system is that it can be found in such states with weight pi, one is 
effectively dealing with a statistical ensemble. Its mathematical description 
is by a density matrix, namely: 

P = ^2Pi\A)(H eM d (C). 

i 

Definition 6 Density Matrix 

Any positive operator p G M^(C) with Tr(p) = 1 describes a mixed state; 
let p = J2j=i r j\ r j)( r j\ be its spectral representation with < r j < 1 and 
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Tj = 1 . Then p defines a positive, linear and normalized functional on 
M d (C) H : 

M d (C) 91^ u p (X) := Tr(pX) = ^M^Wi)- t 3 - 1 ) 

3 

Remark 1 When the density matrix is a projector onto the corresponding 
state vector, which is its eigenvector with eigenvalue one, the state is said to 
be pure, hence the spectral representation of the density matrix reduces in to 
a one- dimensional projection. 

The set of all the density matrices of a quantum system S , is denoted 
by S(S). Since any convex combination of two density matrices is a density 
matrix as well, therefore S(5) is a convex set. Its extreme points, which can 
not be written as a convex combination of other density matrices, correspond 
to pure states. 

Example 4 Block Sphere 

There is a very useful geometrical representation for density matrices p G 
M2(C) of qubits: the Bloch sphere. Taking the two eigenvectors |0) and |1) 
of the Pauli matrix 03 as the axes of this sphere: 



z 




y 



Then any normalized vector in a two-level system can be represented 
in terms of these two vectors: 

= cos0|O) + e <0 sin0|l). 

It corresponds to a point on the surface of this sphere, which can be also 
represented using polar coordination by a new vector ^ 6 tf: 

= (ai, 0:2, CK3) = (sin 29 cos (p, sin 29 sin 0, cos 20) . (3.2) 
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The vector a is known as Bloch vector. In general a mixed density matrix 
p G M 2 (C) can be represented in terms of Bloch vectors and Pauli matrices, 
called Bloch representation: 

1 3 

P = ^2 + ^2^^), (3.3) 

i=0 

where cti = Tr(p<7;). 

The positivity of p is equivalent to asking that ||7^|| 2 = Y2^=i a 1 — 1- Note 
that pure qubit states correspond to Bloch vectors tt with norm one, i.e. to 
the points on the surface of the sphere. Mixed states are instead represented 
by points inside the sphere where \a \ < 1. The maximally mixed state 4^ 
corresponds to the center of the sphere, while orthogonal states, in the sense 
of the Hilbert- Schmidt product, are connected by diameters. 

Therefore a density matrix is a Hermitian, unit-trace operator with non- 
negative eigenvalues. For pure states Tr(p 2 ) = 1, while for mixed states 
Tr(p 2 ) < 1. Given two density matrices, there are different methods to judge 
about their mixedness. 

The mixedness measure on a Hilbert space with dimension d, is defined 
to be [53] : 

M = ^ I (l-Tr(p 2 )). (3.4) 

This quantity is known as linear entropy and is for pure states and maximal 
for maximally mixed states, i.e. cj — -|. 

However, the most reknown entropy is the von Neumann entropy [75] . 

Definition 7 von Neumann Entropy 

Given a density matrix p G 'B(S), the von Neumann entropy is defined: 

d 

S(p) = -Tr(plogp) = -X>logA, (3.5) 

i 

where Aj are the eigenvalues of p and d is the dimension of the system. 

Evidently this quantity is zero only for pure states, as the projectors have 
only eigenvalues equal to or 1. For the maximally mixed states po = 
the von Neumann entropy is maximum: 

s (po) = ~X^ log ^ = S^ logrf = logd - 

i i 

This shows that the von Neumann entropy can measure the mixedness of 
density matrices, and we have < S(p) < log d. 
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3.2 Composite Systems 



In Quantum Information, one very frequently deals with systems consisting 
of several subsystems, called multi-partite systems: S = Si + S 2 + ... + S n . 
The Hilbert space of S is the tensor product of the Hilbert spaces of its 
subsystems: M® n = (^)" =1 Hj. On tensor product Hilbert spaces, operators 
act as tensor product of matrices. 

Definition 8 Matrix Tensor Product 

Consider a bipartite system described by the Hilbert space <C dld2 = C dl <8> 
C d2 , together with the matrix algebra M dl (C) for the first party and M d2 (C) 
for the second one. Thus the algebra of the whole system consists of the 
tensor product of the form flffi [73]/ : 



A®B 



( aiiB cli^B ■■■ ai n B\ 
\0"m,iB a m ^B ■■■ a m ^ n B J 



AeM dl (C),BeM d2 {C). 



When H® n = (g)™ =1 C d \ the algebra of H® n is M® n (£) = (g)™ =1 M di (C). 



The opposite of combining the subsystems to get a bigger one, is partial 
tracing. By tracing over a subsystem Sj with respect to an orthonormal basis 

e one re duces a density matrix p G M® n (C) to another density 
matrix: 

n 

pj := Tr,p = ^(tf'Vl^) e Mf ^(C) = (g)M di (C), (3.6) 

A; i=l 

acting on Cf n " 1} = (g)?=x C*. 

In the following, we will consider only finite-dimensional bipartite systems, 
i.e. S = Si + S 2 , with H 02 = Hi <g> H 2 , Hi = C dl and H 2 = C da . Therefore 
for the whole system we have H = C dl <8> C^ 2 = C dl(i2 . Density matrices are 
positive P12 G M dlrf2 (C) of trace 1. The partial trace yields: 

d 2 

M dl (C) 3 Pl = Tr 2 (p 12 ) := ^(/>i 2 |^> (3.7) 

i=i 
di 

M d2 (C) 9 p 2 = Tn(pi 2 ) := £>J|pi 2 |eJ>, 

i=l 
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where pi (resp. p 2 ) is called reduced density matrix for the system Si (resp. 
for the system S 2 ), the sets {\e\)} d l x and {\ff)}f= 1 are orthonormal basis of 
C dl and C d ' 2 respectively 

The reduced density matrices describe the states that two parties locally 
possess. Their meaning is easy to understand in terms of local observables, 
of the first party, say: A <g> id. The mean value of such a local observable, of 
the first party alone, with respect to the two party state p i2 reads 

Tr 12 (p 12 A® 1) = X>!|A ( X»i 2 |^>) k 1 ) = TtMPi) , 
i=l j=i 

showing that it is completely specified by the reduced density matrix p 1 of 
the considered party. 

For bipartite systems, we have the following lemma: 



Lemma 1 Schmidt decomposition 

Every vector € C dl ® C d ' 2 , can be represented as fj^j -' 

M 

l« = E r J 1) K (1) )®l r f ) ) ( 3 - 8 ) 
i=l 

where M < min{(ii, <i 2 } ; and \\r\ and {1?^ are orthonormal sets 

of vectors in C dl and C d2 respectively. The non-negative scalars {r{\ are such 
thatYZ ir f = l. 



Proof: Without loss of generality let us assume di < d 2 - Given two ONB 

"C I e «) }i=i an d {\fj)}j=ii an y vec tor ipu G C dl <8> C d2 can be written as: 

|« = ^c lJ |e J 1 )®|//), Cij eC. (3.9) 

Its corresponding density matrix is pi 2 = \1p12} ("012 1 • The reduced density 
matrix pi is obtained by partial tracing over the second subsystem, which in 
its spectral representation is: 

Pi = Tr 2 (|0 12 )<0 12 |) = XXV^Krfl. (3.10) 

i=i 



Choosing {\r\ as an ONB for C dl , we can re-write (3.9) as: 
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l« = £X>ik W >®|/;> 
i=i j=i 

= X>< {1) >® 



(3.11) 



where = S/=i c «jl/j)- Note that the vectors \<p\ z> ) are not in general 

orthonormal. Using this new ONB, partial tracing over the second subsystem 
of P12, will be read as: 



(2)\ 



di 



Pi = Tr 2 (|^ 12 )(^ 12 |) = IrfWW^?)- (3-12) 



Comparing this with (3.10), we see that: 

<*FW> = v* (1> - 

After normalization we have: 

l^(2)\ _ « UX)\Jfi\ 



Replacing this in to (3.11), we get: 



i« = Ev r ? ) K (1) )®i r ? ) )- 



i=l 



Which ends the proof. 



(3.13) 



□ 



The coefficients are called Schmidt coefficients, and their multiplicity is 
known as Schmidt rank. Note that using the Schmidt form of a given vector 
^12, the reduced density matrices are: 



M 



Pi 



Tr 2 (p 12 ) = £r l (1) |rJ 1) )(rf ) |, 



M 



P2 = ^{pu) = Y.r? ) VT ) W?\ 



i=l 



one sees that they have the same eigenvalues r- , with the same multiplicity. 



3.3. SEPARABLE AND ENTANGLED STATES 



31 



3.3 Separable and Entangled States 

Consider a vector state \^) G C dl <g> C d2 , which is a tensor product of two 
vector states: |^) = <8> ^2), where G C dl and |*0 2 ) <E C d2 . The 
corresponding density matrix is of the form: 

p* = |*)(*| = IVlXV'll ® |-02> <-02 1 = Pi ®P2, 

where pi )2 are pure states projections. As we see the density matrix corre- 
sponding to a product vector state is itself a tensor product. Using this fact 
we can give a definition for a separable state: 

Definition 9 [30] A bipartite pure state p is separable if and only if it is 
the tensor product of the density matrices of its subsystems: 

p = pi®p 2 . (3.14) 

Otherwise it is an entangled pure state. 

For a pure separable state, tracing over one subsystem gives a projector 
over the other subsystem, i.e. pure states: 

Pi = Tr 2 (pi 2 ) = Tr 2 (pi ® p 2 ). 

Examples of entangled pure states are the Bell states which are maximally 
entangled. 

Example 5 Bell States 

The following four symmetric vectors in C 2 £g> C 2 are known as Bell basis 
vectors LJO^ -' 

= -L(|00)±|11», (3.15) 
= i=(|01>±|10». (3.16) 
For instance consider the density matrix corresponding to the p$+ : 
P $+ = ^(|0)(0|®|0)(0| + |1)(1|®|1)(1|) + 

|(|0><1| <8> |0><1| + |1><0| <S> |1><0|), 
reduced density matrices are: 

1 2 

pi = Tr 2 (p$+) = p 2 = Tr!(p $ +) = — , 
which are maximally mixed states. 
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The following theorem provides a necessary and sufficient condition of 
separability for bipartite pure states, based on their Schmidt decomposition: 

Theorem 2 The pure vector state \1p12) is separable if and only if only one 
Schmidt coefficient is different from zero. 

Therefore the pure states, whose Schmidt decomposition contains more 
than one Schmidt coefficient, or in other words, their Schmidt rank is greater 
than one, are entangled. 

Example 6 According to Schmidt decomposition the following state is sep- 
arable: 

1^12) (lfai\ = (4>1 1 ® 1^2) (V>2 I 

One can easily see that the von Neumann entropy is zero for the density 
matrix p\ 2 , as well as the reduced density matrices p\ and p 2 : 

S(pn) = S( Pl ) = S(p 2 ) = 0. 



Example 7 The Schmidt rank of Bell states in (3.15) is 2, and their Schmidt 
coefficients are 

The symmetric states in (3.32) are another example of entangled states, 
which Schmidt rank d. 



3.3.1 Shannon entropy and von Neumann entropy 

Example [5] is a clear indication that entanglement embodies correlations that 
have no classical counterpart. Indeed, the von Neumann entropy is for quan- 
tum states what the Shannon entropy [19] is for classical states, that is for 
probability distributions of stochastic variables. Let X = {xj}^ =1 be a d- 
valued stochastic variable whose outcomes Xj occur with probabilities Pi, 
then the Shannon entropy of X is defined by 

d 

H{X) = -^PdogPi . (3.17) 

i=l 

One thus sees that the von Neumann entropy of density matrix p is noth- 
ing else than the Shannon entropy of its spectrum. However, the Shannon 
entropy of a composite system, described say by two stochastic variables 
X\ = {x^}^ and X 2 = {x^}^^ is always larger than the Shannon en- 
tropy of each of its two parties, that is 



H(X!, X 2 ) > max{FpTi), H(X 2 )} , (3.18) 



3.3. SEPARABLE AND ENTANGLED STATES 



33 



where H(Xi,X 2 ) is the Shannon entropies of the joint probability distri- 
bution 7Ti )2 = {p(xf\x^)} and if(X lj2 ) are the Shannon entropy of the 
reduced probability distributions 7Ti = {p(^^)}f=i, ^2 = {p(^f ^j^li with 



j=l i=l 
This can be seen as follows: let us introduce the conditional probability 



'i \ x j > 



that the outcome for X\ be x\ given that the outcome Xj has been obtained 
and the conditional entropy of X\ given X 2 , namely the following convex 
combination of Shannon entropies 

H{X l \X 2 ) = -f^pixf^f^pix^lxf^logpix^lxf) (3.21) 

j=l i=l 

= H(X 1 ,X 2 )-H(X 2 ) . (3.22) 
Then, the result follows because H(X 1 \X 2 ) and H(X 2 \X 1 ) are > and 

H(X 1 ,X 2 ) = H(X 2 ) + H{X l \X 2 ) = H{X 1 ) + ifp^^) . (3.23) 



Remark 2 Unlike the Shannon entropy, the von Neumann entropy of sub- 
systems can be larger than that of the total system. This fact indeed signals 
entanglement but, as we shall see, is not always implied by entanglement. 

One could now guess that quantum states either carry entanglement or 
classical correlations. It is not so: there can be separable states which are 
not simply classical correlated. 

In order to understand this fact one has to introduce the so-called Quan- 
tum Discord of a given state. One starts by recalling that, using the condi- 
tional entropy, one can define the mutual information: 

: X 2 ) = if(Xx) - if (Xi|X 2 ) (3.24) 

which measures the amount of information that is common to both classical 
stochastic variables. It can also be rewritten as: 



3(X 1 : X 2 ) = H(X 1 ) + H(X 2 ) - H(X U X 2 ). (3.25) 
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One can extend the notion of mutual information to quantum systems 



; To do so, starting from (3.25), one can simply replace the probability 



distributions of individual systems by reduced density matrices of subsys- 
tems: 

3a-.b(pab) = S(p A ) + S(p B ) - S(pab)- (3.26) 



However this generalization, in the case of expression (3.24) for classical mu- 



tual information, is not trivial; in fact, the definition of conditional entropy 
of quantum systems involves quantum measurement processes. Expressing 
the state of system A, when the state of system B is known, is equivalent 
to having done a set of measurements on system B. Therefore the post- 
measurement state of system A, after performing measurement Hf will be: 



PA\nf 



1 

— ( 

Pi 



X 



ufpABid ® nf ), Pi = Tr( PB nf ). 



iB 



(3.27) 



Given a complete set of orthogonal measurements II = {Ui,i = 1, ...,d}, 

:IT = l d , ] 

defined as [501 



^jilj = Id, UiUj = SijUi, on system B, the conditional entropy will be 



S(p A \U B ) :=^p,S(p A |nf). 



(3.28) 



With this definition of quantum conditional entropy, Z in (3.24) reads: 



3a-.b(p 



AB 



S(p A ) - S(p A \Il 



B\ 



(3.29) 



Unlike in the classical case where the two expressions for mutual infor- 
mation, 3 and 5, are equal, in quantum case, the two expressions are not in 
general the same [501 EH]- Their difference depends on the set of measure- 
ments performed on system B which can be removed by minimizing it over 
all possible sets {II B }, and define the minimum as the Quantum Discord of 
the given state pab'- 

Da-.b(pab) = mm[3 a-.b(p - AB) a-.b(pab)} 
{ri s } 

= mm[S(p B ) + S(p A \n B ) - S(pab)] ■ (3.30) 
{ri s } 

The Quantum Discord measures the non-classical correlations whhc are present 
in a quantum state: it is for separable states with only classical correla- 
tions, and maximal for entangled states which are extremely non-classical 
correlated. However, there exist separable, that is non-entangled, states with 
non-zero Quantum Discord. 
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3.3.2 Mixed entangled states 

For mixed states the notion of separability extends as follows: 

Definition 10 A mixed state p is separable if and only if it cannot be written 
as a convex combination of product states: 



P 



E 



PiPl ® P2> 



E 



Pi 



1, Pi > o. 



(3.31) 



Otherwise, it is an entangled mixed state. 

Example 8 Consider a bipartite system S = S± + S2, where the Hilbert 
spaces of the subsystems are of the same dimension. Consider the symmetric 
states 



*i> = ^X>->®b>> 



(3.32) 



where \j), j = 1,2, ...,d is any fixed orthonormal basis in C d . The corre- 
sponding projectors are: 



(3.33) 



^/so m i/izs case reduced density matrices are maximally mixed: 

1 d 

pi = p2 = 



i=i 



Id 
d 



which means that \^ d + ) is entangled. We shall refer to this set of symmetric 
states in the future. 



Example 9 Werner States 

Werner states in C 2 £g> C 2 are constructed as follows I5b\ IZfij : 



Pc 



1 — a 



4 



-l0 ® lr 



< a < 1. 



(3.34) 



Where is one of the Bell states, introduced in (3.15). In the basis where 
(T3 is diagonal, Werner states can be represented by the following matrix: 



Pa 



1 

4 



(\-a \ 

1 + cv -2a 

-2a l + a 

\ 1-aJ 
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The above range for a is to guarantee the positivity of p a . Werner states are 
entangled for a > |. The reduced density matrices are: 

Oa)l = {Pa)2 = ^ 

As we see, also in this case the reduced density matrices are maximally mixed 
and therefore their von Neumann entropy is maximum. Indeed, by checking 
how the function A a = S(p a ) — log 2 behaves: 




one might then think that entangled Werner states are detected by checking 
whether their entropy is smaller than the entropy of one of the parties. We 
shall later see a counterexample to such a claim. 

However the situation in case of pure states is completely clear thanks 
to Schmidt decomposition, deciding if a given mixed state is entangled or 
separable is still an open problem, even in case of bipartite systems. In the 
following chapters we will study how to partially answer this question. 




Chapter 4 



Positive and Completely 
Positive Maps 



Physical transformations of Quantum systems are described by linear maps, 
either on the space of states or, dually, on the algebra of observables. If a 
linear map on the convex space of states describes a physical transformation 
(that, of course, transforms physical states into physical states), then it must 
transform density matrices into density matrices. It must thus be positivity 
preserving and trace preserving. However, though necessary, positivity is not 
sufficient to guarantee the physical consistency of a linear map. Via coupling 
the system affected by the transformation with a suitable unaffected ancilla, 
the existence of entanglement between the two subsystem requires the map to 
be also completely positive. In this section we will see how these constraints 
can take part in the challenge of detecting entangled states. 



4.1 Positivity and Complete Positivity 

Definition 11 Positive Map A linear map A : M m (C) — > M n (C) is posi- 
tive if and only if M„(C) 3 A[X] > for all X e M m (C) J22 [57]/. The set 
of positive maps is denoted by V. 

Example 10 Trace Map 

The Trace Map, introduced in (2.7), is a positive map, as it is the sum 



of the eigenvalues of an operator. Therefore if M^(C) 3 X = ^2 i Xi\Xi)(xi\ ) 
Xi > 0, then 

Tr : X i—). Tr(X)l d = (J^)^ > 0, Vie M d (C). 
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Example 11 Transposition Map 

Transposition map is another example of positive maps: 

T d : M d (C) 9lMl T 6 M d (C) (4.1) 

with respect to a fixed orthonormal basis. Given a matrix X = [Xij], its 
transposed is X T = [Xji\. Since transposing an operator does not affect its 
spectrum, it preserves positivity, which means that it is a positive map. 

Example 12 Transposition Map on M 2 (C) 

Using the Pauli matrices as an orthonormal basis set in M 2 (C), similarly 
to the Trace Map in Example [3[ we can also write the Transposition Map 
as ffl: 

3 

M 2 (C) 9lM T 2 [X] := ^e a a a Xa a (4.2) 

a=0 

where e a = 1 when cr ^ 2, whereas e 2 = —1. In terms of the map S a as in 
(2.12), the Transposition Map is of the following form: 

3 

T[X] = ^2e a S a [X]. (4.3) 

One can thus check that, in the standard representation, T[a a ] = a a , 
a 7^ 2 and T[a 2 ] = — o 2 , which is exactly what transposition does in the 
chosen representation. 

Definition 12 K-positivity 

A linear map A : M m (C) — > M n (C) is k-positive if and only if id^ £g> A : 
M fc (C) <S> M m (C) -> Af fc (C) ® M n (C) is a positive map, where k e N. V k 
denotes the set of k-positive maps f75[ \6h 



Definition 13 Complete Positivity 

If a linear map A is k-positive for every k E N, it is called completely 
positive. The set of completely positive maps is denoted by CP . 

A linear map A is a CPU map if it is both CP and unital, that is A[l] = 1. 

The following theorem fully characterizes the completely positive maps: 

Theorem 3 Kraus representation 

A map A : M m (C) — > M n (C) is completely positive if and only if there 
exists a set of operators {Ki}, such that J^3| \H\ \3jj: 

A[X] = K\XK U VX e M m (C). (4.4) 

i 

The operators M nxm (C) 3 Ki : C™ — > C m , are called Kraus operators. Fur- 
thermore, such a map is CPU if and only if K^ = 1. 
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Theorem [3] concerns the action of maps on operators (Heisenberg repre- 
sentation); however, one is often interested in how states transform (Schrodinger 
representation). This fact can easily be derived by duality. 

Definition 14 Dual map 

For any map A : M m (C) -»■ M n (C), its dual map A T : M n (C) -»■ M m (C) 
is defined as FJ]/; 

Tr(A[X]p) = Tr(XA T [p]) (4.5) 
for every X e M m (C) and p G §„. 

The Kraus representation of the dual map A T : M n (C) — > M m (C) will 

be: 

i 

Therefore, one can see that CPU maps and trace-preserving maps are dual: 
A CPU <J=^ A T Trace -preserving. 
Consider a CP map A : M n (C) — > M n (C), with its Kraus representation 



as in g6}. By choosing a Hilbert-Schmidt ONB set {i^}™ =1 G M n (C), the 



Kraus operators -K^ can be written as: 



Therefore the Kraus representation of A will be: 



i j,fc=l 
n 2 

i,fc=i 

Where fjk = J2i c ijCik are the entries of a positive matrix F = [fjk] £ 



M n 2(C). The representation (4.7) of a CP map A is also known as Canon- 



ical Kraus Representation [7j. Indeed the Kraus form (4.6), is obtained by 
diagonalizing this general form. 

One should also note that the set of Kraus operators {Ki} is not unique, 
however, each two sets of Kraus operators corresponding to the same CP 
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map, can be related by a unitary matrix. Indeed, given a canonical Kraus 
form, of a CP map A: 

n 2 

A[X] = djkFjXFl, < G = [g jk ] e MAC). (4.8) 
j,k=i 

Using the diagonal form of G in terms of its eigenvectors: 

n 2 

G = ^2gi\gi){gi\, gi>0, 
i=i 

the matrix elements are: 

n 2 

9jk = {j\G\k} = ^gigijM- 
i=i 



Replacing this in (4.8), we get: 



A[X] = 9imWkF 3 XFl 
j,k,i=i 

9 9 9 

1=1 j=l k=l 



Any hermiticity and trace-preserving linear map A : M^(C) — > M^(C) can 
be written as [27] : 

d 

A|>] = ogltfpF, , (4.9) 

where Ylij=i a ijFi Fj = This can be shown as follows: consider the d 2 
elementary linear maps X — > Tij \X\ = Fi X Fj and define on the linear space 
£ of maps A : M^(C) — > Md(C) the scalar product 

« A 1 |A 2 »=rf 2 Tr(id,®A 1 [P^]tid (i ®A 2 [P^) , (4.10) 



where P+ is the completely symmetric projection (3.33). It turns out that 
<< Fij\Fki »= SikSjf, therefore, being orthogonal, the Tij constitute an 
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orthonormal basis in C. As we shall soon see, the scalar product defined 
above is strictly related to the notion of Choi matrix. 

When the map A is completely positive the matrix A = [a^] is positive, 
while in case of positive maps, the matrix A = [ay] is no longer positive, 
but with separating the positive and negative eigenvalues, one can write a 
positive map A as the difference of two completely positive maps: 

am = *i k Ip k i - E \ a i\ K }p R i- ( 4 - n ) 

a ; >0 a;<0 

Based on the complete characterization of the structure of CP maps, the 
following concept provides a simple test to determine if a given map is CP 
or not [52] . 

Definition 15 Choi-Jamiolkowski isomorphism 

Given a linear map A : M n (C) h- > M m (C), its Choi matrix is the matrix 
on C n <g> C m defined by [M EI 

M nxm (C) = M n (C) <g> M m (C) 3 C A = id„ x A[P| n) ] , (4.12) 

where projects onto the completely symmetric state 

1 - 

|* (n) ) = - Y h> ® \ei) , {|e 4 )}™ =1 ONB in C n . 
i=i 

It is useful to see how the Choi matrix of the map A looks like: 



/A[|l)(l|] A[|l)(2|] ■■■ A[|1)H]\ 
A[|2)(l|] A[|2}(2|] ■■■ A[|2)(m|] 

VA[|m)(l|] A[|m)(2|] ■■■ A[\m){m\)) 



The Choi matrix, C\, is uniquely defined by A since its entries C&iij, P q) 
with respect to an orthonormal basis of the form | e^) <g> \ fj) G C™ <S> C m are 

C HHm) = ( e * ® fi\ Ca K ® /,) = -(/i|A[| ei )(e p ]|/ g ) . (4.13) 

It turns out that the map J : A i— > J [A] = Ca from the vector space C n ^ m of 
linear maps onto the matrix algebra M nxm (C) is invertible map and defines 
the so called Jamiolkowski isomorphism. In fact, any A : M n (C) i— > M m (C) 
is determined by its action on matrix units \ei)(ej\ which is in turn com- 
pletely characterized by the matrix elements in \J^.13 ). Therefore, given 
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C G M nxm (C) this uniquely determines a linear map Ac '■ M n (C) i-> M m (C) 

(/,-|A c [|e i >(e p |]/,> = C(y, pq) . 

Equipped with these tools one can readily prove the following 

Theorem 4 Every linear map A : M n (C) i— > M n (C) positive semi- 
definite Choi matrix C\> is completely positive and vice versa J7g[ [^/. 

Proof: By definition of complete positivity, id n <E> A is positive if A is com- 
pletely positive, and thus Ca = id n (g) A[P( n )] > 0. 

Vice versa, if M n 2(C) 3 Ca > 0, then it can be spectralized as 

n 2 

The eigenvectors \cj) G C n £g> C n can be conveniently rewritten using the 
completely symmetric vector as follows: 

p,q=l p=l 

with the matrix Cj G M n (C) defined by 

(e q \ Cj \e p ) = nc?j 9 . 

Therefore, one can rewrite 

n 2 

C K = Y; 1 n® K 3 P + )l n®K} , K 3 = ^Lj CjEMniC), 
J'=l 

Then, by the Jamiolkowski isomorphism, A[X] = 52?=i -^j AT can be writ- 
ten in the Kraus-Stinespring form and is thus completely positive. □ 

Remark 3 In general, positive maps A : B(M) i— > -B(H) on bounded opera- 
tors are completely positive if idd <g> A /or all d > 2. However, in the finite 
dimensional cases, when B(M) = M n (C), the above result shows that one 
needs check the positivity of id n <8> A, only. 
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Example 13 Consider the Trace Map $ : M d (C) i — > M d (C): 

$(X) = ^Tr(X)l d . (4.14) 

It is easy to see that this map is CPU; indeed: 

$[ld] = Id- 
It's Choi matrix is: 

= ^]£|i>01® ^(1001)1001 

TTie positivity of Choi matrix shows the Trace Map is completely positive. 



Definition 16 Block Positive Operator 

Given a matrix A G Mj(C) ; using any fixed orthonormal basis, we can 
write it as A = Ylij ®Aj, where Aij G M d (C), for VI < i,j < d. Then 
matrix A is said to be block-positive if and only if for all \<p) and G C d , 
the following holds: 

(0®VW®V>> = (0|(^K#)|z)(j|)|0> 

({MMM ^\A 12 \^) ■■■ (i>\A ld \^)\ 
(MAM ■■■ MA 2 M 



> 0. 



In another word, block-positivity of A is equivalent to positivity of 
for G C d . 

Theorem 5 A linear map A : M d (C) — > M d (C) is positive if and only if its 
Choi matrix is block-positive /ffPj 

(ip®<p\l d ® A[P d + ]|^®v3) > 



for every \ip) and \(p) in C d . 
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Proof: Given a positive map A : M^(C) — > M^(C), by definition, it pre- 
serves positivity of matrices in M^(C): 

(0A[|^)(^|]|0)>O VhA>,|0> GC d . (4.15) 

Given an ONB {\i)}f =1 e C d , then any vector \ip) is = X^iV^K), where 



ipi = (i\ip). Using this (4.15) is 



= (^®0ix;ii>oi®A[|i)O'i]iV'*®0) 

= (^®0|l d ®A[P+]^*®0) 

= (^®0|C A |^*®0) 
> 0. 

Which shows that the positivity of A is equivalent to the block positivity of 
its Choi matrix and vice versa. □ 



Example 14 Consider the Transposition Map T : M<f(C) — > Md(C), in (4.1). 
It's Choi matrix is: 

ci = lJ2\i)(j\®\j)(i\. 

Which coincides with the definition of the Flip operator: 

d 

V =H l*>01 ® b'X*l = T)P+, (4.16) 

which is also that 

Since V 2 = idd, its eigenvalues are ±1. 

For any vector \<p) and G C d , we have: 

(^j ® ip\c T \ij ® ip) = (ip®<p\(%2\i){j\ ® b'XW®^) 
= X^<^|i)O-|0><Vli)<*l^> 
= £<#X#>EWXV>b'} 

« i 

= I(^I0)| 2 
> 0. 
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Since Ct has negative eigenvalues, Ct 0, Transposition Map is not com- 
pletely positive. 

Example 15 Reduction Map 

In the next chapter we will introduce a number of positive maps, which 
are known to be able to detect some entangled states, one of them is called 
Reduction Map [32], which has the following form: 

A R = Ti d - id d (4.17) 

The Reduction map is positive: as the trace is the sum of the eigenvalues, 
which are non-negative as p is a density matrix, therefore Tr^fp] > 0. Let \ipi) 
be eigenvectors of p with eigenvalues \, then using the spectral decomposition 
of p, we have: 

3 « 

Since ^ \ > Xj, Ar[p] is a positive matrix. On the other hand, the Choi 
matrix of this map is: 

C Ar = id d <g> A R [P+] = l d 2 - dP+ 
1 S 

where = - \i) (j I ® |«) 0' I • 

Thus, C\ R has one negative eigenvalue 1 — d and is not positive semi-definite. 
Therefore, the Reduction Map is not completely positive. 

In what follows, we shall explain the results given in [US], which is exten- 
sively used in the last chapter of this thesis. 

Theorem 6 Let Tr be the Trace Map. Then, given any positive map <ft : 
M d (C) Mrf(C) can be written as 

0=^(Tr-0 cp ), (4.18) 

where c > is a positive constant and (f> cp : M d (C) —> M d (C) a completely 
positive map. 
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Proof: 

Let be a self-adjoint map, that is (/)[A*] = (<f>[A]y, then its Choi matrix 
is self- adjoint as well. Therefore can be written as difference of two 

positive matrices, i.e. = — where C^C^ = 0. Let fi > be the 

smallest positive number such that 

iii>c* = c}-c;. 

The above inequality is satisfied by choosing c = ||Ct||, which is the greatest 
eigenvalue of C£. 

Now consider the following positive matrix: 

C£ = l--^>0. (4.19) 

The completely positive map <j) cp of which is the Choi matrix reads 

<Pc P = Tr - -<j>, (4.20) 
fx 



Now, re-arranging the (4.20), we have: 

0=^(Tr-0 cp ). (4.21) 
Note that the only assumption was to consider the map to be self-adjoint. 



Therefore any self-adjoint map can be written as in (4.21). Taking into ac- 
count that the set of positive maps is a subset of self-adjoint maps completes 
the proof. □ 

Example 16 We have seen in Example^ that the transposition on M 2 (C) 
can be represented as 

1 3 

M 2 (C) 914 T[X\ = - ^ s a S a [X], (4.22) 

o=0 

S a [X] = cr a Xa a , e a = (1, 1, -1, 1) . 
Therefore, on M 4 (C), 

1 3 

M 4 (C) 914 T[X] = ^J2 £ « £ ? S ^ X ^ ( 4 - 23 ) 

a,/3=0 

S a p[X] = a a pX cr a p , a a p = a a ® . 
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Notice that the transposition is not written in the Kraus form (4-6) as the 
products e a ep = — 1 whenever a ^ (3 = 2 or (3 ^ a = 2. 

On the other hand, the completely positive trace map Tr on M 2 (C) has 
the following Kraus representation: 



1 3 

Tr[X] = -Y,S a [X] , XeM 2 (C) , 



(4.24) 



a=0 



Then, on M^iC), the trace map has the Kraus form (4-6): 

1 3 

M 4 (C) 914 Tr[X] = - S aP [X\ . 

a,B=0 



(4.25) 



We would like to write the transposition map as in the form suggested by 
Theorem 

A = /i(Tr — A CP ). (4.26) 



When A is the transposition map on M 4 (C), the CP maps <pcp i n (4-^6) 
are easily found: the linear map 



r a, /3=0 



E a Ep 



a/3 



is of the form (4-6), thus CP, if and only if fx > 1. The smallest choice, 
\l = 1, yields the completely positive map 



A 



CP 



2 ( Yl Sa2 + Yl s ^ 

a^2 /3^2 



(4.27) 



4.2 Entanglement and Positive Maps 

As we have seen the structure of completely positive maps is wholly un- 
derstood via the Kraus representation. Moreover Choi's theorem provides 
a technique to sort them out. Instead, as soon as one is dealing with 
positive maps, the situation becomes much less clear due to the lack of a 
complete structural representation of them. As we will see, only in low 
dimensions, when d\ x d 2 < 6, the structure of positive maps is under con- 
trol [75, 63J. While in higher dimensions these maps are only partially un- 
derstood [201 EH ESI EDI [21 SB]. 
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From a physical point of view, it turns out that a linear map which de- 
scribes a physical transformation needs not to be only positive but completely 
positive since the eigenvalues of a system in statistical interpretation 
represent the probabilities, when a linear map describes a physical trans- 
formation, it has to preserve the positivity of the spectrum of the physical 
system, otherwise a negative eigenvalue representing a probability brings up a 



contradiction. As we have seen, by Definition 13 any extension of completely 
positive maps does this job. However this might not be true in general if the 
transformation is described by a positive, not CP, map. While the extension 
of positive maps preserves the positivity of separable states, it can fail when 
the subsystems are entangled. Consider a bipartite separable state which is 
a convex combination of tensor product of states its subsystems: 

P = Y,PiP 1 i®P 2 i ^>°> I> = L ( 4 - 28 ) 

i=l i 

The action of id& <g> A when A is a positive map gives: 

p' = (id fc g> A)[p) = Y^PiiMpl] ® A[p']- 

i=i 

The positivity of A guarantees that A(pf) and therefore p', remains a positive 



state. But if the state p is not separable it cannot be written as (4.28), and 
consequently, for a generic p and k, (id& <g> A)[p] is positive only if A is 
completely positive. 



However, though positive maps cannot be used to describe physical trans- 
formations, they provide a new powerful theoretical approach to detect the 
entanglement. Indeed, if a given state p does not remain positive under the 
action of an extended positive maps, it must be entangled: 

(idfc <S> A) [p] ^ =>■ p entangled. 

The prototype of positive, but not completely positive map is transposi- 
tion. In the following example we show that the extension of this map fails 
to be k-positive already when k is equal to 2, which means that it is not a 
complete positive map. In other words, though a positive matrix remains 
positive after being transposed, it might not be positive after being partially 
transposed. 

Example 17 Take one of the maximally entangled Bell states, introduced 
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in 



(3.15), |$ + ) = ^=(|00) + 1 11)), its density matrix is: 





/ 1 








1 \ 


1 














2 
















V i 








1 / 



The action o/1.2<8>T on this state is to transposing each block of the matrix, 
the result is the following matrix: 



(l 2 ®T)p 





/ 1 








o\ 


1 








1 





2 





1 










\o 








1 / 



which has a negative eigenvalue — |. 

According to Partial Transposition, the states can be divided into two 
classes: those which remain positive under this map, called PPT states, and 
those which turn out to have at least one negative eigenvalue after being 
partially transposed, known as NPT states. 

The following class of positive maps has been used to classify these maps 
in low dimensional Hilbert spaces. 

Definition 17 Decomposable Map 

A positive map A : M d (C) — > M d (C) is decomposable if and only if it 
can be written as [6$ : 

A = ri + T 2 o T d (4.29) 
where Fx and T 2 are complete positive maps. 

Theorem 7 All the positive maps A : M dl (C) — > M d2 (C) where d\ x d 2 < 6 
are decomposable |75| [ftffi . 



Example 18 The reduction map Ar : M d (C) — > M d (C) in (4.17) is decom- 
posable J3EJ . 

The Choi matrix of this map reads: 



1 



dPj 



C Ar = 1®A r [P+] 
1 d2 

where P d + = - ^ N) 0' I ® I*) 01 ■ 
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Its partial transposition with respect to the second subsystem yields: 



(l d ®T d )o(l d ®A R )[P+] 
(l d ®ToA R )[P+] = ^(l d *-V). 



(4.30) 



Since the flip operator V, introduced in (4.16), has eigenvalues ±1, it im- 
plies that C^ 2 R has no negative eigenvalue and therefore is a positive matrix, 
whence the corresponding map A' = TAr is completely positive, according to 
Theorem [^} Furthermore, since Ar = TA' where A' is CP, results decompos- 
able. 



In the next chapters we will see how and under which conditions positive 
maps can detect the entanglement. 



Chapter 5 

Detecting Entanglement 



In the previous chapter we have seen formal necessary and sufficient condi- 
tions for separability, practically, establishing whether a given state is sepa- 
rable or not is in general a very hard task. Indeed apart from PPT criterion 
in low dimensions, the available theorems do not provide constructive tools 
to detect entangled states in general. In the present chapter we will review 
some of the entanglement detecting methods. Though none of them is able 
to detect all entangled states, they nevertheless work for particular classes of 
states. 

5.1 PPT Criterion 

In the last chapter, we have introduced Transposition Map. Since this map is 
not CP, it can be used to detect entangled states. In fact, partial transposi- 
tion provides one of the most easily applicable tools to detect entanglement. 
Unfortunately in higher dimensions, being PPT is no longer a sufficient con- 
dition for separability. 

Theorem 8 Horodecki Criterion 

A bipartite state p G SdxcZ is separable if and only if for any positive map 
A : M d (C) Mrf(C), (id d <8> A)p is positive fifTf. 

As this must hold for any positive map A, so does it for Transposition 
Map: 

Theorem 9 Peres Criterion 

If a bipartite state p G Sdxd is separable, then (idd <8> T)p > 
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Obviously when (id^ ® A)p ^ 0, the state is NPT, and according to the 
above theorem, is entangled. As we have already mentioned, Transposition 
Map is an exhaustive entanglement witness in low dimensions: 

Corollary 1 PPT Criterion A state p acting on C 2 <g>C 2 , C 2 <g>C 3 or C 3 <g> 
C 2 is separable if and only if its partial transposition is a positive matrix f31^ 



Proof: 

If p is separable, Theorem [8] tells us that its partial transposition p T2 is 
also positive. 

The sufficiency is proved as follows: In Theorem [7] we have seen that all 
positive maps A : M^C) — > Md 2 (C) when d\ x <i 2 < 6 are decomposable 
and can be written as: 

A = r 1 + r 2 oT, 

where are CP maps. The extension of this map acting on p gives: 

(id d ® A)[p\ = (id,, ® ro[p] + (idd ® r 2 )[ P T2 ]. 



The first term is positive because I\ is a CP map and by definition (idd ® 
> 0. If p is PPT, then p T2 > 0, then the second term remains positive 
as well, due to complete positivity of T 2 , which completes the proof. □ 



Example 19 All Bell states introduced in (3.15), have negative partial trans- 
position and are NPT entangled. 



Example 20 In (3.34), we have introduced the Werner States in C 2 £g> C 2 : 

1 — a 



4 



<a<l. 



Their partially transposed density matrix is: 



T 2 



P, 



Its eigenvalues are: 





(I 


— a 








-2a \ 


1 







1 + a 








4 










1 + a 







V 


-2a 








1-aJ 



Xi — X 2 — A3 



1 + a 



-, A4 



1 — 3a 
4^ ' 
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As we see the first three eigenvalues are positive, while the fourth one when 
| < a < 1 is negative. Therefore: 

- 1 1 

— - < a < - =^ p a is separable, 

o o 

- < a < 1 =>• p a is entangled. 

In Example^ we have seen the behavior of A a = S(p a ) — log 2, where we 
have mentioned, though it might look like that von Neumann entropy might be 
able to detect the entangled Werner states, however, it only detects entangled 
states for 0.75 < a < 1, but it fails to detect the entangled Werner states for 
l<a< 0.75. 

5.2 Entanglement Witness 

The question of whether a given density matrix p of a higher system is sepa- 
rable or not is very difficult to answer due to the difficulties of reconstructing 
the whole density matrix in terms of tensor product of its subsystems, i.e. to 
write p = pi | ipj ) | <8> | ipf ) ( ipj | • However the fact that a set of such separa- 
ble states is convex and closed, makes it possible to apply the Hahn-Banach 
separation theorem. This opens a new approach to detect entanglement, 
called entanglement witnessing. 

Theorem 10 Hahn-Banach Theorem 

Let X be a finite- dimensional Banach space and S a convex compact set 
in it. Let xq be a point in X such that xq ^ S. Then there exists a hyperplane 
which separates x from S FTB} 0/. 
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Let the set S se „ correspond to the set of separable states which by defini- 



tion is convex, compact and closed. Therefore the Hahn Banach theorem 10 



says that there must exists a hyperplane separating this set from an entan- 
gled state which lies outside. This separator is called Entanglement Witness, 
which is a Hermitian operator J3TJ EI]- 



Definition 18 Entanglement Witness 

A Hermitian operator W acting on H = C dl <g> C d2 is called an entangle- 
ment witness if 



1. 3p G S dlXd2 such that Tx(pW) < 0. 

2. (<p ® ip\W\<p ® ip) > V|v?) e C d \ \tp) e C d2 . 

3. Tr(W) = 1 



The second property shows the block-positivity of the non-positive oper- 
ator W, while the last one asks only for normalization. There are immediate 
consequences: 



Theorem 11 p G §d X d is entangled if and only if there exists an entangle- 
ment witness such that Ti(pW) < JffI] /. 



Theorem 12 a G S^xd is separable if and only ifTr(o~W) > for all entan- 
glement witnesses W WT^ . 



Definition 19 Decomposable Entanglement Witness 

Let P and Q be positive operators acting on C d <g) C d , then W is a decom- 
posable entanglement witness if and only if |76|/ 

3a, b > such that W = aP + bQ Tl 



where Ti is Partial Transposition. 



Theorem 13 An entanglement witness, EW, is non-decomposable if and 
only if it detects PPT entangled states J^ffi . 



Theorem 14 p G SWd is PPT entangled if and only if there exists a non- 
decomposable entanglement witness such that Ti(pW) < J31\ \7b^ . 
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Given two entanglement witnesses W\ and W 2} the latter one is said to 
be finer [IB] if all the entangled states detected by W\ are also detected by 
W 2 . An EW, Wopt is optimal [121 EE] if there exists a separable state p 
such that Tr(pW op t) = 0. It means that the hyperplane corresponding to 
W opt is tangent plane to E> sep (for illustrations see Fig (5.2)). Since optimal 
entanglement witnesses detect more entangled states than those non-optimal 
ones, they are highly interesting. Indeed the set of separable states can be 
completely characterized by the set of optimal EW [TT] . 




It has been shown [16] that a non-optimal EW can be optimized, but 
one should note that none of these entanglement witnesses can detect all the 
entangled states. 

Both positive maps and entanglement witnesses provide necessary and 
sufficient conditions for separability. Indeed to some extent their behavior is 
quite similar toward entangled states; there exists a positive map and an EW 
for each entangled state. This brings up the conjecture that these connections 
can be explained through an isomorphism. 

Definition 20 Jamiolkowski Isomorphism 

Jamiolkowski isomorphism says that there is a one-to-one relation between 
an EW as an operator and a positive map A : M d (C) — >■ M d (C) [39]: 

A w < — >W A e M d2 (C) 

Given an operator W\ acting on M d 2 (C) and choosing an orthonormal prod- 
uct basis set {|e, <S> fi)} in H = C d ® C d , the corresponding linear map 
A w : M d {C) — ► M d (C) is defined LJ^: 



A w {p] = Tr 2 {Wp T >}. 
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Conversely from a given map A, the operator W\ is obtained: 

W A = (l d ®A)\^)(^\ 



where \ip^} is defined in (3.32) 



The following properties hold [HJ |39l [76] 

• W > if and only if Aw is completely positive. 

• W is EW, if and only if A w is positive. 

• W is decomposable if and only if A w is decomposable. 

• W is non-decomposable if and only if Aw is non-decomposable. 

5.3 Reduction Map 

In the previous chapter, we have introduced a positive (not CP) map : 
M d (C) — > M d (C), called Reduction map: 

A R [p] = Ti d [p]-p. (5.1) 

Below we will see how this map can be used to detect the entanglement: 

Theorem 15 Reduction Criterion 

A bipartite state p\i acting on H = C dl £g>C d2 with dixd 2 < 6, is separable 
if and only if JM/ 

idd <E> P2 - P12 > 0. 

For states acting on higher dimensional Hilbert spaces, the criterion is only 
necessary for separability. 

This map however can not detect PPT entangled states, as we have shown 
in the last chapter, that this map is decomposable. 

5.4 Range Criterion 

The range criterion is another approach to detect entanglement. This method 
is based on finding the product vectors in the range of the states. The range 
of p is the set of all the vectors for which there is another vector \<p) such 
that: 

\iP)= P \4>). 



5.4. RANGE CRITERION 
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Definition 21 Product Vector 

Let H = (QTLi C* fre an m-partite Hilbert space. A vector G H is 
called product vector if and only if it has the form \ip) = \4>i) where 

|&)eC* 

Note that a product vector is not necessarily separable with respect to 
all the partitions of Hilbert space. For example consider the vector = 
^|001) + |111). This is a product vector when we look at it a bipartite state, 
i.e. in C 4 <S> C 2 it can be written as \ip) = ^(|00) + |11)) <g) However in 
three-partite system in C 2 <8> C 2 <E> C 2 the vector is no longer a product 
vector. 

Theorem 16 Range Criterion 

If a bipartite density matrix p i2 acting on Hilbert space H = C dl Cg>C d2 , is 
separable then there exits a set of product vectors {e^ <S> fi} where | e^) G C dl , 
\fi) G C d2 , where |/*) is complex conjugate , such that {ei ® fi} spans the 
range of p i2 and {e^ ® /*} spans the range of WM- 

Example 21 The following PPT state in C 3 Cg> C 3 has been shown jSffi to 
be entangled: 
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l+a 
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where < a < 1. 

Note that the range criterion is neither stronger nor weaker than PPT 
criterion, as it can detect PPT entangled states, while there are NPT states 
which do not violate this criterion either. 

There are some states which extremely violate this criterion, such states 
are called Edge states [^5] 153]. 

Definition 22 Edge state 

An edge state 5 is a PPT entangled state such that, for all product vectors 
|e<g>/) (|e) G C dl and \ f) G C d2 ) and e > 0, 5 - e\e g) f)(e <g> f\ is not positive 
or does not have a positive partial transposition [34\ l3R 
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Obviously the edge states lie between PPT entangled states and NPT 
states. 

Example 22 Indeed the first PPT entangled state in C 2 eg) C 4 found in [3^] 
is an example of an edge state: 
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Vl-6 2 
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1+6 
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where < b < 1. 



Using an edge state S, one can construct a non-decomposable entangle- 
ment witness as follows |16] : 

Let us denote the range and kernel of 5 by R(5) and K{5) respectively 
Let P and Q denote the projector onto K(5) and K(S T2 ) respectively and 
define 



where a 



a(P + Q T2 ) 
1 

Tr(P + Q) 



(5.2) 



Define also 



inf (e®f\W 5 \e®f). 

|e®/> 



(5.3) 



Theorem 17 fjm Given and edge PPT entangled state 5, then W oc Ws— el 
is a non-decomposable entanglement witness, where e and Ws are defined 
m Q and 



Proof: 

To show that W is an entanglement witness, we have to show that it 



satisfies the requirements in (18). The block positivity of W requires that 
for any G C dl and |0) G C^ 2 we have: 



(tp ® (f)\W\ip ® <f)) > 
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The first two terms cannot vanish simultaneously, as 5 is an edge state: If 
{ip ® 4>\P\tp ® (p) = then ® (f>) G R(S), then by definition of an edge 
state, \ip®4>*) R(5 T2 ), which implies that \ip®(jf) G K(5 T2 ) and therefore 
(tp <S> 4>*\Q\ip £g> 4>*) = 1, and vice versa. Therefore (ip <g> <p\Ws\ip <S> 4>) > for 
all |V) and \<f>). 

On the other hand W detects 5: 



Tr(W5) = a(Tr(P5) + Tr(Qj)5)) - eTr(15) 



The first two terms are zero due to fact that P5 = Q5 T2 = 0. Since W de- 
tects a PPT entangled states, by (13) it is a non-decomposable entanglement 
witness. □ 



5.5 Unextendible Product Basis 

Searching for the product states in the range of a density matrix was leading 
into a new method called Unextendible product basis, mentioned for the first 
time in [5]. 



Definition 23 Unextendible Product Basis 

Consider a multipartite quantum system H = (^)^ x C di with m parties. 
An orthogonal product basis (PB) is a set S of pure orthogonal product states 
spanning a subspace M s of H. An uncompletable product basis (UCPB) is 
a PB whose complementary subspace H -1 i.e. the subspace in H spanned by 
vectors that are orthogonal to all the vectors in M s , contains fewer mutually 
orthogonal product states than its dimension. An unextendible product basis 
(UPB) is an uncompletable product basis for which M s contains no product 
state. 



As an example of product basis which does not span the whole space, one 
can take |0) <8> |0) and |0) ® |1). As we see these two orthogonal vectors form 
a PB which spans a two dimensional subspace in C 2 ® C 2 . 



Example 23 Consider the following example of UPB given in JEB$ , known as 
Tiles. It is a set of five product vectors in C 3 ®C 3 , spanning a 5-dimensional 
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subspace: 



lV>o) 
lV>i) 
1^2) 
1^3) 

|V>4> 



|0>®^(|0>-|1» 

|2>®^(|1>-|2» 
1 



y/2 
1 

1 



(|0>-|1»®|2> 
(|1>-|2»®|0> 



(5.4) 



(|0> + |1> + |2»®(|0> + |1> + |2». 



One sees that as all of these vectors are mutually orthogonal and each three 
of them span a full 3-dimensional Hilbert space, therefore it is not possible to 
find another product vector which would be orthogonal to all of them. 

For the bipartite Hilbert space H = C d £g> C d , when d is even, one can 
construct UPB as follows [21]; Let's denote the set of d orthonormal vectors 
by |0), \d — 1), then the set of vertical tiles are constructed as follows: 



d/2-1 



\lpmn) 



\n) ® \u m>n+x ) = \n) <g> ^2 ujm \i + n + 1 mod ( 5 - 5 ) 



d 



m 



j=0 



— 1 and n = 0,...,d — l, 



where to = e* 47r//d . Similarly, the horizontal tiles are defined to be: 

d/2-1 



v) 



m 



\wm,n) ® \ n ) — yj w* m \j + n mod d) <8> |n), 

j=0 



(5.6) 



1,.., 



1 and n = d — 1. 



The following theorem illustrates the role of the UPB in the challenge of 
detecting entanglement. 

Theorem 18 The state that corresponds to the uniform mixture on the space 
complementary of a UPB {ipi : % = 1, n} in a Hilbert space of total dimen- 
sion d 



P 



1 

— (i-EfovWil) 



d — n 



(5.7) 



is a PPT entangled state JMj. 
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Proof: 

Partial Transposition maps a set of UPB into another set of UPB, while 
leaving the identity unchanged, this shows that the state in (5.7) is indeed a 
PPT state. On the other hand, p by construction contains no product state 
which means that it is entangled, all together it is a PPT entangled state. □ 



5.6 Realignment 

We first introduce some new notions which will be used in the following. 

Definition 24 Reshaping 

Consider an m x n matrix A = [a^] where are its entries. By putting 
the elements of each column one after another, one obtains a vector, which 
can be associated to this matrix. Conversely, a vector of length mn can be 
turned into an m x n matrix PR 



A 



m,n 



«2,1 «2,2 



0>l,n\ 
«2,n 



< > ~ct 



an, 



5 a m X, 



(5.* 



7 



Let X be an mn x mn block matrix and a^i = (an, • • • , a n i, • • • , a nn ) T 
the vector corresponding to the block kl of size n. The realigned matrix X R 
of size m 2 x n 2 , is defined by: 



X 



exampl 


e let us 


take 


m = 


n = 2: 


X = 


( ^1,1 

^2,1 


Xl,2 
X 2 ,2 


Xl,3 
X2,3 


Xl,4 \ 
%2,4 


^3,1 
\ ^4,1 


^3,2 
^4,2 


^3,3 
^4,3 


x 3,4 
^4,4 / 



The vector corresponding to the first block is an 
and the realigned matrix is: 



(5.9) 



(xn, X21, X12, X22) 
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X 



R 



%2,1 %l,2 x 2,2\ 

#3,1 x 4,l x 3,2 x 4,2 

^1,3 ^2,3 x l,4 X 2,A 

^3,3 ^4,3 ^3,4 ^4,4/ 



In [7J , reshaping a matrix is performed by aligning its rows one by one in 
a vector, instead of its columns: 



A, 



/ aii 01,2 

^2,1 0-2,2 



&2,n 



<; S> a — (an, «12, " " " , Oln; • • • , 0"mn) 



Now any mn x mn square matrix can be reshuffled by reshaping each row 
into a square matrix and putting them block after block into a new matrix. 
The new reshuffled matrix, X R obtained from matrix X in (5.9) is: 



R' 



( 


X\2 


x 2,l 


X 2,2 \ 




X X,A 


^2,3 


X2,4 


^3,1 


^3,2 


X± t \ 


^4,2 


\ X 3,3 


2^3,4 


^4,3 


^4,4 / 



X 



In general, given a matrix X G M^(C) ®M^(C) = M^C), let us consider 
the orthonormal basis consisting of tensor products \e m <g) f^) G C d2 , where 
{l e -m)}m=i an d {|/^)}f=i are orthonormal bases in C d . Further, let us denote 
as in [7J by 

X^ = (e m ® / M |X|e n ® /„) . (5.10) 



Then, the matrix elements of reshuffled matrix X R are given by 



vR' 

j\ mf_i 



(5.11) 



Remark 4 Note that the result of both methods are equivalent up to a per- 
mutation. In both cases after reshaping a matrix into a vector, the Hilbert- 
Schmidt product /?/ will be read as scalar product between two vectors: 

(A\ B) = TxA^B = (a'\b'). 



Reshuffling the tensor product of two matrices turns to be very simple. 
Consider two matrices X G M P (C) and Y G M q (C). We denote the vectors 
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assigned to each matrix by vec(X) and vec{Y). Now consider their tensor 
product: 



(x X \ Xl2 
\x p i X p 2 

( xwVn 
£112/21 

X\\Vq\ 



X2p 



X 



I'!' 



yu yi2 

2/21 V22 



V2q 



\y q i y q 2 ■■■ y qq J 

XllVlq \ 
XllVlq 



X12Y 



xipY 



XllVqq 



X2lY 



X22Y 



x 2p Y 



(5.12) 



\ Xp\Y Xp2Y • • • XppY ) 

By the above method the vector corresponding to the first block is: 
at = (xuVu, • • • , x n y lq , • • • , x u y qq ) T = x u Y, 



where is the vector corresponding to the matrix Y (see (5.8)). Putting all 
these vectors as rows of a new matrix, we have the reshuffled matrix of this 
tensor product: 



R 



fxn(^) T \ 

x lp (^y 

[x PP (^) T ) 



We now show that reshuffling a density matrix can be useful to detect its 
entanglement. 

Theorem 19 Reshuffling Criterion If a bipartite density matrix p 12 is 
separable, then the trace norm of its reshuffled matrix will not be increased fTQ 
Ef, i.e.: 

1 1 Pl2 1 1 Tr < ||/?12 || Tr = 1- 

Proof: 

If P12 is separable, then it can be written as pyi = J2iPiPi ® Pi where 
J2iPi = 1- Consider the pure separable state: pu = p 1 <S> p 2 - As we have 
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seen in the last remark, by reshuffling this state we have: pf 2 = /?[/?] T . But 
in this case, p 1 and p 2 are projectors and therefore their norm is equal to the 
norm of the corresponding vectors. We have: 



llpf 2 || < ll?[?] T ll < \\pi\\ = Hull W = i- 

By properties of the norm, we can see that: 



l/4ll<J>ll(Pi®^)*ll<Z>t<i- 



□ 



This theorem provides a necessary condition for separability, therefore 
if one finds that for a given density matrix p, after reshuffling ||p R || > 1, 
this means that the state is entangled. Note that p R is not Hermitian. It 
was proved [601 EE] that realignment is independent of Partial Transposition, 
therefore none of them is neither stronger nor weaker than the other one. 



Example 25 As an example one can consider the density matrix constructed 



as in theorem (5.7), using the UPB given in example (5.4), called tiles: 



i=0 



However we know that this state is entangled but we can also see it by direct 
calculation of the norm of \\p R \\ = 1.32. Since \\p\\ < \\p R \\ it reveals the 
entanglement of this state using the reshuffling method. 



Chapter 6 
cr-diagonal States 



In the following chapter, we consider bipartite states that are diagonal in 
the basis generated by the action of tensor products of the form 1 2 »® ap, 
cfft = ®™ =1 cr w , on the totally symmetric state |\&+) G C 2 ™ <g> C 2 ™. We 
first characterize the structure of positive maps detecting the entangled ones 
among them; then, the result will be illustrated by examining some entangle- 
ment witnesses for the case n — 2. Further, we will show how, for pairs of two 
qubits, being separable, entangled and bound entangled are state properties 
related to the geometric patterns of subsets of a 16 point square lattice. 



6.1 cr-diagonal and lattice states 

In the following we shall freely call entanglement witness both positive maps 
A such that id^ <S> A[p en t] ^ and their Choi matrices M\. 

We shall consider a bipartite system consisting of two parties in turn 
comprising n qubits; the corresponding matrix algebra M 2 2n(C) is linearly 
spanned by 4™ tensor products of the form a ft := £8>™ =1 o" Mi = cr m £g>cr M2 £g>- ■ • a^ n . 

Given the totally symmetric state, already introduced in (3.32): 

i*+> = ^X>')®b> 

with d = 2 n , the vectors 

|^> := 1 2 „ ® a p \^) G C 2 " ® C 2 " , (6.1) 
form orthogonal projectors 

P fl := I^X^I = (1 2 » ® ^)|*+")<*f |(la« ® ^) • (6.2) 
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Orthonormality follows since 

t=l i=l 

Definition 25 T/ie c/ass o/ bipartite states we shall study will consist of 
states of the form 

P = J2 r * P t ' < ^ < 1 , £ t-0 = 1 , (6.3) 

i/iai zs diagonal with respect to the chosen orthonormal basis in C 2 ™ £g> C 2 ™: 
we s/iaiZ caiZ i/iem a-diagonal states. 

A particular sub-class of states of two pairs of two qubits, n = 2 and 
°~pi = °~a® o~p = o~ a/ 3 G M 4 (C), are considered in jH \5\ E] and called lattice 
states. In this case the orthonormal basis vectors and their corresponding 
projectors are the following: 

|*^> = l 4 ®a a/i |* + }eC 16 , 

P a p = \^afi){^afi\ = (U ® CT Q/ 3)P 4 + (1 4 ® (Ta/3), 



where Pf is defined in (3.33) 



For them partial results concerning their entanglement properties have 
been obtained. We shall tackle them again in the following: the class of 
lattice-states is defined as follows. 

Definition 26 Lattice States 

Taking a subset I C L 16 of cardinality Nj, then the corresponding lattice 
state pi is defined by: 

pi = jf E p <# ■ ( 6 - 4 ) 

1 <*,Pei 

Let Li 6 denote the set of pairs (a, 0), where a and (3 run from to 3: it 
corresponds to a 4 x 4 square lattice, whereas the subsets 

C a := {(a,/3) e Li 6 : = 0, 1,2, 3} 

Rp := {(a J) eL X6 : a = 0,1,2,3} (6.5) 

correspond to the columns and rows of the lattice, respectively. 



If compared with those in (6.3), the lattice states are uniformly dis- 



tributed over, and thus completely characterized by chosen subsets of Li 6 . 
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Example 26 As an example consider the lattice states associated with the 
following subsets: 



' (-Poo + Pq2 + Pll + Pl3 + P20 + P2I + P23 + P< 



32 1 



Ni = 8 



N T = A : 
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X 
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X 






X 


1 




X 


X 







X 




X 









1 


2 


3 


+ P11- 


- P22 + P33) 


3 








X 


2 






X 




1 




X 









X 













1 


2 


3 



(6.6) 



(6.7) 



In the following, we shall be interested in investigating the entanglement 
properties of such states. As an appetizer, we consider whether the von 
Neumann entropy of any pj can be smaller than that of its reuced density 
matrices. Indeed, we have seen that, if so, pi would be entangled. By tracing 
over the first or the second party we get as reduced density matrices 



Tr 2 ,ip/ 



with maximal von Neumann entropy S(p} 



I4 
4 

(l),(2h 



2 log 2. On the other, 



end each lattice state is already spectralized with equal non-zero eigenvalues 
1/iVj; therefore, S(pi) = log Nj > 2 log 2 whenever Nj > 4. 



Furthermore, if we consider the Reshuffling Criterion in Theorem 19 



using the relations (|5.10|) and (|5.11|), we obtain as reshuffled lattice states 

(pr 



J TThpu 

nu 



nm 



1 

AN~! 



^ (n\<r a p\v) {n\a a p\m) 



whence 



Pf 



1 

AN~! 



E 

(a,/3)e/ 



a a(3 ® a al3 , 



where we chose the ortonormal basis vectors in ( |5.10 ) to be the ones used 
in writing the Lattice States, namely |ej) = \fi) = Thus, using the 



68 CHAPTER 6. a-DIAGONAL STATES 



properties of the trace norm, one gets 

whence for Nj > 4 no entangled lattice state can be detected by the chosen 
criterion. 

We now proceed to examine the lattice states from the point of view of 
partial transposition. Positivity under partial transposition (PPT-ness) of 
lattice states pj is completely characterized by the geometry of I [I] . 

Proposition 1 A necessary and sufficient condition for a lattice state pi to 
be PPT is that for every (a, (3) G Lie the number of points on C a and Rp 
belonging to I and different from (a, (3) be not greater than iVj/2. In terms 
of the characteristic functions Xi{ a -,f3) = 1 if («,/?) G I, = otherwise, a 
lattice state pi is PPT if and only if for all (a, (3) G L 16 : 

3 3 pj- 



Example 27 Consider the following lattice states: 
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3 
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3 



(6.8) 



(6.9) 



These states do not remain positive under partial transposition and are there- 
fore entangled: in the first state p$, the row and column passing through 
the point (2,2) ^ / contains 4 > 5/2 points in I, while in the second one 
3 > 4/2 = 2 points. 
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Example 28 By the same criterion, the following two states are instead 
PPT, 



Nt = 6 : 
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x 3 
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X 






x 2 


X 
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x iV 7 = 8 : 1 






X 


X 
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3 



(6.10) 



They need not be separable as in lower dimension; indeed, a sufficient crite- 
rion devised in [3] show them to be entangled. 

Proposition 2 A sufficient condition for a PPT lattice state pi to be en- 
tangled is that there exists at least a pair (a,/3) G L 16 not belonging to I such 
that only one point on C a and Rp belongs to I. Equivalently, pj is entangled 
if there exists a pair (a, 0) G L w , (a, (3) ^ I , such that 



0=5^/3 0=5^a 



In both patterns of the states in Example (28), it is the point (0,0) ^ / 
which satisfies the sufficient criterion. 

Example 29 Unfortunately, this criterion fails in the case of the PPT lattice 
state characterized by the following subset 



Nj = 10 
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3 



(6.11) 



is 



Indeed, the only candidate point to fulfill the criterion in Proposition ^ 
(2,2); however, it belongs to I. 

Luckily, a refined criterion [6] based on [12] shows it to be entangled. 
Proposition 3 A PPT lattice state pi is entangled if the quantity 
k i U = X/(a,f©2)+ X/(/i©2,/3), 

where © denotes summation modulo A, is such that kj U = 1 for a column 
and a row R u ®2, independently of whether (/x©2, v®2) belongs to I or 

not. 
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Example 30 Thus, the state in (29) fulfills the sufficient condition k®° = 1. 

Another example of a lattice state detected by the last proposition is the 
following: 



11 : 
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As it is apparent from these examples the various criteria for entanglement 
are closely related to the structure of the subsets that define the lattice states; 
in the following we will try to clarify this correspondence. 



6.2 Entanglement detection for cr-diagonal states 

In this section, we will show that positive maps that be entanglement wit- 
nesses for states of the form (6.3) can be sought in a particular subclass of 
them. 

Any linear map on A : M<f(C) i— > M^(C) can be written as [3]: 

d 2 -i 

M d (C) 3X^ A[X] = £ K GiXG) , 

i,j=0 

where the matrices Gi G M^(C) form an orthonormal basis with respect to 
Hilbert-Schmidt scalar product, namely Tt(G\ Gj) = 5ij and the coefficients 
Xii are complex numbers. 



In the present case, the normalized tensor products J— constitute such a 

V2 n 

basis in M 2 n(C), whence linear maps A : M 2 n(C) h-> M 2 n(C) can be expressed 
as 

M 2 n (C)3 X ^ A[X] = J2 V SplX] , S m [X\ = apt X a, . (6. 12) 



ere 



The next one is a simple observation based on (6.1) and (6.2). 

Lemma 2 A o-diagonal state p = rp, is entangled if and only if th 
exists a positive map A as in (6.12) such that ^jiii r fi < 0- 

Proof: According to Horodecki Criterion in ([8]), the state p = rp Pp 
is entangled if and only if 

Tr (id 2 n ® A[Pf ] p) < , 



6.2. ENTANGLEMENT DETECTION FOR a-DIAGONAL STATES 71 



for some positive map A : M 2 «(C) i— y M 2 n(C). 

Using the orthogonality of the vectors \^/p) in (6.1), p is entangled if and 
only if there is a positive map A such that: 

Tr (id 2 n <g> A[Pf ] p) = V (^u\p\^it) = X ^ r t < ' 

p,,T) fi 

□ 

The lemma indicates that the class of diagonal positive maps of the form 
A = J2ft -V Spii might suffice to witness the entanglement of states of the 
form p = J2ft r fi Pp.- What is needed to show that it is indeed so is the 
following result. 

Lemma 3 Given a positive map of the form A = - , the diagonal 

map A diag = ]T\ is also positive. 

Proof: From Choi's theorem ([5]), we know that a given map is positive 
if and only if its Choi matrix is block positive. Considering the map A = 
J2ft v -V* 7 Sp$, the positivity assumption is equivalent to 

Y,^vWM)^\°M) > V|^),|^)eC 2n . (6.13) 

Given a pair \i[)),\<f) G C 2 ™, consider another pair \ipSi) = &5i\i>), Iv 9 ^) = 
(Tg^ip), where cr^ denotes the tensor product where pj = for j ^ i and 
A*i = o~i 7^ 0. Inserting the new pair into (6.13), we get 



a^m x 



o Vj ® (o h o Vi a^ (g) <T Vj \ip) > . (6.14) 



jX,v j=l j=i+l 

i-1 

j=l j=i+l 



Consider pi ^ uf, because of the Pauli algebraic relations, one can always 
choose in such a way that 

o- Si (j^(J5. i = a H and o 5i u Vi u &i = -a Vi , 



whence all the terms in (6.14) corresponding to the chosen pair of indices 
(pi, Vi) contribute with 

i—l n i—1 n 

j=l j=i+l j=l j=i+l 
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Then, adding inequalities (16.13) and (6.14) yields 



By applying the very same argument for pairs (pi, Vi) with varying index i, 
one cancels all contributions from p^ 7^ i/j and remains with 

^ A# H^Wk/ilv) > V |^), |^> G C 2 " . (6.15) 

This, by Choi's theorem (|5| amounts to the positivity of the diagonalized 
map A diag = A/z/2 S^. □ 

The previous result allows us to conclude with 

Proposition 4 Entangled p = r^P^ can be witnessed by diagonal posi- 
tive maps A = Y Jf i A/i/j Sffi. 

Proof: By the previous lemma, diagonalizing a positive map 

A = A^i? Sfj 



always yields a positive map Adiag = Yla Spp. Then, from Lemma ^ it 
follows that either the entanglement of p is witnessed by an already diagonal 
map or, if by a non-diagonal one, also by the map obtained by diagonalizing 
the latter. □ 



Remark 5 From the above Proposition we only get that entangled states of 
the form (6.3) can be witnessed by diagonal maps Adi a g — 5-^5 A/z ^fijii the 
main problem is of course to characterize the coefficients in such a way 



that (6.15) be satisfied and thus A d i ag be positive. 



Expression (4.25) for Trace map in M 4 (C) can be extended to the trace 



operation on M 2 n(C), 



Tr 



In the second part of chapter 4, we have seen that positive maps can be related 



to the complete positive maps. This relation was given in Theorem (4.26), 
which allows us to write a positive map as: 



A = MT-A cp ). 
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Therefore diagonal positive maps A diag are also related to completely positive 
diagonal maps 

P 

by setting 

A^ag = fJ- £ - A m) % ' > ' ( 6 - 16 ) 

P 

and asking for block positivity of its Choi matrix, which is 

MWIMV»)>o, V|^),|^)eC 2 ", 

(^l(Tr-ASr a )[|^)^|]|^) = l-X]^|(^|a^)| 2 > V \tp), G C 2 " 

(6.17) 

Using these observations, we get necessary and sufficient conditions for 
the separability of a-diagonal states. 

Proposition 5 A o -diagonal state p = r p Pp. is separable if and only if 
for all sets of 4™ positive real numbers Xp > 0, such that 

XixMaMl 2 <1 V|^),|^)GC 2 " , (6.18) 

P 

it holds that 

p 

Otherwise, if a set of 4 n positive real numbers > can be found that 
satisfy (6.18) and such that 

£ *prp > ^ > ( 6 ' 2 °) 

p 

then the o -diagonal state p = Y2p r pPp entangled. 

Before tackling the case of lattice states, that is of a-diagonal states with 
n = 2, as a simple application, we consider the case of only two qubits, n = 1, 
for which we know PPT-ness to coincide with separability. 

Example 31 In the case n = 2, a-diagonal states have the form 

3 3 

p = £vP /t , r M > , £r> = 1 

H=0 fj,=0 
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and the Pa's project onto the Bell states, already introduced in (3.15): 



,T v , T 2v |oo) + |ii) 



V2 



l*i> 



|00)-|11) 



1*3) 



|Q1) + 1 1Q> 

01) - |10> 
y/2 ' 



(6.21) 



Under transposition T[a^\ = er M; = (1, 1, —1, 1) and under partial trans- 
position the projection Pq goes into the flip operator V\ip ® if) — \ip <g> ip): 



1 1 

id®T[P } = -V=-J2 VpP*, ^ = (1,1,1,-1) . 



Therefore, the action of partial transposition on a o-diagonal state yields 



1 3 

id ®T[p] = - JJr^lS^Vl®^ 

fi=0 

1 3 

= 2 ^ ^ r M P M 

^3 3 

1 3 
1 3 



(6.22) 



a=0 



Definition 27 In the above expression, the following construction has been 
employed: given (a,fi), a,// = 0,1,2,3, is i/ie unique index from to 

3 suc/i i/iai o~ a a^ = rjail &[a,/j], where rja/1 is a phase ±1 or ±i, they can be 
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considered as non-zero elements of the following matrices: 



/O 1 0\ 

10 

i 

\0 -i 0/ 

/O 1\ 

i 
-i 

\1 0/ 



/O 1 \ 

-z 

10 

\0 i / 

/l 0\ 

10 
10 

\0 1/ 



Because of the Pauli algebraic relations, the symbol 
properties that can be used to derive (6.22): 

[a, ji] = [/!, a] , 



enjoys the following 



[a, lA = 7 [a, 7] = P => \M, 7] = « ■ 
27ms, [a, ■] is a one to one map /rom £/ie set (0, 1, 2, 3) onto itse//. 

.Since positivity under partial transposition identifies all separable states 
of two qubits, a a-diagonal state is separable if and only if r^ < - for all 
u = 0,l,2,3. 

On the other hand we have seen that all positive maps acting on C 2 ® C 2 
are decomposable, therefore using (4.26), ( |4.22[ ) and (4.24), we can write any 
diagonal positive map A diag : M 2 (C) >-)■ M 2 (C) as 



diag 



J2^ } S a + J2^S a oT 

3 1 3 

E A ^^ + 2E ep^SaoSp 

a=0 a,/3=0 

=0 

3 1 



7=0 



0=0 



A 7 j o 7 , 



7=0 



where Xa' 2 ^ are, according to Theorem (4.4), positive numbers. Then, the 
coefficients 
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can always be made positive and thus A = X!a=o ^ a ^ a completely positive, 



by choosing p large enough. Then, they fulfill the condition (6.18) that cor 



responds to A$ og being positive. We can thus consider when and whether 



inequality (6.19) which reads 

7=0 ^ 7=0 \ 13=0 / 

is satisfied. By choosing Ay = for all 7 and Ay = 5 ia , one gets 
3 

^ r[ a ,p] ep = r [afi] + r [ai i] - r [a>2 ] + r [afj3 ] = 1 - 2r [Q>2 ] > 

/3=0 

and £/ius, fry varying a, r M < - /or a// = 0, 1,2,3. Vice versa, if r M < - 
/or a// /i = 0, 1, 2, 3, one obtains 

33 33 

7 =0 /3=0 a=0 (5=0 

3 



3 

A « 2) (l-2r M ) >0. 



a=0 



6.3 Lattice states 

In this section we shall restrict ourselves to the lattice states pi G M 16 (C) 



in (6.4), namely to uniformly distributed a-diagonal states with n = 2 



Proposition (J5J now reads 

Corollary 2 A lattice state pi is separable if and only if 

E ^ ^ t (6 - 23) 

(«,/3)e/ 

/or a// choices of 16 coefficients A a p > swc/j £/ia£ 

E A^|(^|o- a/3 |^)| 2 < 1 V |V) e C 4 . (6.24) 

(<*,/3)6/ 
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Otherwise, if a choice of positive coefficients exists that satisfy (6.24) and for 
which 

, Ni , 

E A ^ > t ' (6 - 25) 

then a lattice-state pi is entangled. 

Before drawing concrete conclusions from this result, we examine the 
entanglement criteria in Propositions ^ and ^ in the light of the diagonal 
structure of witnessing maps which is the main result of the previous section. 



Example 32 The states considered in Example (28) were found to be entan- 
gled by showing that they do not remain positive under the action of id ® T t 
where 



r* = 000 (*) 5*00 + ^ (goi(t)S(H + gn(t)SjD 



(6.26) 



i=i 



with 



9oo(t) 
9oi(t) 

9io(t) 



1 + 3e" 4t 3 + e 



-At 



1 + 3e 



-At 



-At 



-It 



3 + e 



-At 



4 4 

was proved to be a positive map from M 4 (C) into itself. This map, expressed 



by means of the notation of the last Example in chapter 4 (4.25), is already 



in diagonal form; the diagonal completely positive maps associated to it by 



St0rmer Theorem (4.26) have the form 

r* 

A cp (t) 



Tr- 

1 

4 

3 

£ 

i=l 



9oo(t)\ c 
I >->oo 

fi ) 

1 _ 9oi(t) 
4 fi 



(6.27) 



9io{t) 



S h 



if) 



with /i which has to be adjusted taking into account that 



1 1 

9oi(t) < — , gait) < — , goo(t) < 1 . 
16 16 
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Then, complete positivity of the map A cp is guaranteed by n>\ which yields 



Aoo(t) = 


1 

4 


9oo(t) 


> o 


X 0i (t) = 


1 

4 


9oi(t) 

n 


>o, 


\o(t) = 


1 

4 


9io(t) 

n 


> o 


\ij = 


1 

4 







These coefficients surely satisfy condition (6.24) as the latter just reflects 
the positivity of the originating map T f ; they also satisfy condition (6.25). 
Indeed, 



M*) 

a ,pei 



4 



(6.28) 



- ^2 [9oa(t))6 a> o6p t o + g O p(t))5 a}0 + g a o(t))8p, ) 
iV 7 1 



(1 - 4t)5 Q ,o<W + t{8a,o£p + <W (6.29) 



a,pei 



For both subsets in (28), the second term in (6.29) is negative due to £2 = — 1. 
Thus, Y^ Q o eI X a p(t) > ^ for small times. 



Example 33 Let us now consider the lattice state in Example (29): 



Nj = 10 : 



3 


X 






X 


2 




X 


X 




1 


X 


X 




X 





X 


X 




X 







1 


2 


3 



(6.30) 



In JM/, it has been shown to be entangled by using the following positive map 

M 4 (C) 3X^$ V [X}= Tr[X]-(T[X\+V[X]) , V[X] = V ] XV , (6.31) 

consisting of the trace map to which one subtracts the transposition map and 
a completely positive map V constructed with a 4 x 4 matrix V such that, in 
the standard representation, 



a+2 P+2 



-<" T . E(i 



V a 2\ + \V2c 



1 . (6.32) 



a^2 
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In this way, 

= 1 - I^WI - ^ MM ^ = i - i 3 - Q , 

where P and Q are orthogonal one-dimensional projections and thus ensure 
the positivity of the map. 

Because ofV, the map $y is non-diagonal in the maps S a p: in order to 
illustrate the content of Proposition (J5j) we diagonalize it. Using 



1 3 

M 4 (C) 914 T[X\ = ^J2 £ « £ P S ^ ' = a <# X a ^ ■ ( 6 - 33 ) 



a,/3=0 

one gets 



The mean value of the Choi matrix of $^ a9 with respect to the lattice state 
in (6.30) reads 



Tr( Pl id®^ a9 [P^ 



i n 



iV/ V2 



v 12\ 



and becomes negative choosing | 1 2 > 1/2 hence revealing the entanglement 

of Pi- 
Proposition associates to $^ a9 completely positive maps of the form 

(b dia9 / 1 1 In l 2 \ 

whose coefficients can always be made positive by choosing \i > 2. The sum 
of the coefficients corresponding to the subset I of the lattice state in (6.30) 
yields 

/or i/ie choice of V\ 2 which exposes the entanglement of pj, in agreement 
with (f6725l) . 
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6.4 Separable Lattice States 

In absence of a positive map that witness the entanglement of a lattice state 
pi, one can check its separability by trying to express it as a convex combi- 
nation by other lattice states that are known to be separable. For some pi 
this is feasible as shown in [6]; for instance, consider the lattice state 



Pi 



\{Pn + Pn + Put + P21 + P23 + P31 + P32 + P33 



According to Proposition Q, it is PPT and is also separable; indeed, the 
defining subset / splits as follows 



3 




X 


X 


x 3 




X 


X 


3 










2 




X 




x 2 








2 




X 




X 


1 




X 


X 


x = 1 




X 


X 


+ 1 




X 




X 




































1 


2 


3 





1 


2 


3 





1 


2 


3 

> 


1 


3 




h 


X 


x 3 




h 

X 




X 


2 








2 




X 




X 


+ 1 






X 


x + 1 



































1 


2 


3 





1 


2 


3 

> 



The 4-point subsets /j are not disjoint, but all points contribute exactly twice 
to /, thence one rewrites 



Pi 



2=1 



in terms of rank-4 lattice states corresponding to the subsets ij. The re- 
sult follows since the criterion of Proposition Q ensures that they are all 
PPT [35]. 



6.4.1 Special Quadruples 

A more general sufficient condition for the separability of lattice states can 
be derived by introducing the notion of special quadruples. 

Definition 28 Special Quadruples 
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A special quadruple Q is any subset of the square lattice L 16 consisting of 
4 points (a, /3) such that there exist \ip), \ip) G C 4 for which 



1 . 



(6.34) 



(o,/3)6Q 



Given a lattice point (a, (3) G L 16; we will denote by Q a p G Q any special 
quadruple containing (a, (5), by Q a/3 the set of such quadruples and by n aj3 
its cardinality. 

Example 34 Consider the lattice state: 

1 



Pa 



{Pq2 + Po3 + Pl2 + Pi 



13 J 



JVj = 4 : 



5?/ choosing \(p) = | H — ) and 



3 


X 


X 






2 


X 


X 






1 


























1 


2 


3 



+ +), where cr 2 |±) = =r|=f) an d 



cr 3 |±) = |=f) ; we see that the set {a 2, 003, oyi, o"i3} satisfies (6.34) 



Example 35 Some more examples of special quadruples are the following 
lattice states: 



3 








X 


3 


X 








2 






X 


2 




X 






1 




X 




, 1 






X 







X 















X 







1 


2 


3 







1 


2 


3 


3 






X 




3 




X 






2 








X 


2 






X 




1 




X 




1 








X 





X 









X 













1 


2 


3 





1 


2 


3 



Since |(v9|cr a/ 3|V')| 2 = 1 if and only if a a p\il>) = r]\(p) with r\ a pure phase, 
a set of 4 points {(<x,-, /3j ) }|=o C ^ is a special quadruple if and only there 
exist lib), \d>) G C 4 such that 



Vj = 0,1,2,3. 



(6.35) 
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Let us focus upon Q 00 , the set of all special quadruples {(0,0), 
(a 2 , /3 2 ), («3, Z^)} containing the point (0, 0). Each of them is obtained from 
the fact that, using (6.35) with (a (3 ) = (00), 



\<p) = Wi 



J = 1,2,3. 



Therefore, the cr a/ g of a special quadruple must commute. Indeed, 



(7/3(75 — 7 7 (7 Q ® (75(7/3 



(6.36) 
(6.37) 



where 77^ are the coefficients ±1 and ±i, introduced in Definition (|27j), such 



that 7q,7/3 = r^aOp and the 16 coefficients e a7 = ±1 are given by the Pauli 
matrix commutation relations, which can be considered as the elements of 
the following matrix: 



(6.38) 



/l 1 1 


1 \ 


1 1 -1 


-1 


1 -1 1 


-1 


v - 1 - 1 





It thus follows that 



:/ 3 , a 7 5 ytp) = if and only if e a ^eps = 1, whence 

Qoo consists of the following 15 special quadruples (omitting the point (0, 0) 
common to all of them) 

{(0,1); (1,0); (1,1)} {(0,2); (2,0); (2, 2)} {(0, 3); (3, 0); (3, 3)} 
{(0,1); (2,1); (2,0)} {(0, 2); (1, 2); (1, 0)} {(0, 3); (1, 3); (1, 0)} 
{(0,1); (3,1); (3,0)} {(0, 2); (3, 2); (3, 0)} {(0, 3); (2, 3); (2, 0)} 



{(1,1); (2, 2); (3, 3)} {(1, 2); (2, 3); (3, 1)} (6.39) 
{(1,1); (2, 3); (3, 2)} {(1, 3); (2, 2); (3, 1)} 
{(1,2); (2,1); (3, 3)} {(1, 3); (2, 1); (3, 2)}. 

Knowledge of <2oo is sufficient to derive the form of all Q a /3] in order 
to prove this fact, let us introduce the following family of maps indexed by 
lattice points (a, (3) G Li 6 : 



r a/ 3 : L 16 h> L 16 , Tapl^S)] = ([«, 7 ], [(3,6}) , (6.40) 
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where the map (0,7) H- [0,7] has been introduced in Definition (27). It 
follows that the maps r a p are invertible: T a p o T a p[(j,5)] = (7, 8). Given a 
subset / = {(ai,fli)} C Li 6 , r aj g[i] will denote the subset {T aj g[(ai, A)]}- 

Lemma 4 T/ie map 9^4 T Q/ g[Q] € Qoo one-to-one. 

Proof: If Q = {(%, 0)}? o is a special quadruple for (a, 0) = (a ,fl ), 
then a a jfa\ijj) = e* Xi |0) for some \ip), \<f>) G C 4 . Right multiplication by a a p 
yields: 



e lx '>o- a p\<f)). 



Then, {([a, otj), [0 Pj]}j=o = T a( g[Q] is a special quadruple for (0,0) exposed 
by the vectors and o~ a p\(j)). The one-to-one correspondence follows from 
the invertibility of the maps r a p. □ 



It is easy to check that all Q G <2oo in (6.40) give rise to lattice states pQ 
that satisfy the criterion in Proposition (|lj) for being PPT; this is also true 
for lattice states corresponding to Q G Q a p: indeed, they are obtained from 
the previous ones by the local action a a p <g> \pqa a p ® 1. 



Lemma 5 Properties of the special quadruples Q 



00 



1. Given two commuting a a p and <r 7< 5 there is a unique a^ u commuting 
with both of them. 

2. Two quadruples cannot have more than one pair (a, 0) in common; 
otherwise they would coincide. 

3. Each a a p commutes with 6 a 1 s different from itself and from ctqo 

4. Each a a p belongs to three different quadruples. 

5. The pairs (7, 5) 7^ (a, 0) belonging to the 3 quadruples with (a, 0) in 
common correspond to anti- commuting a 's. 



Proof: The first and second property follow directly from (6.37, they im- 
plies the third and fourth ones, while the last one is a consequence of the 
fourth property. □ 

The following Lemma shows that they are indeed separable as dictated 
by the general result in [35]. 
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Lemma 6 Let I C L 16 ; all positive coefficients \ a p satisfying inequality (6.23) 
are such that Y^( a p)eQ -"W — ^ / or a ^ special quadruples Q C J . 



Proof: Recall the trace map (4.25): 

1 3 

M 4 (C) 3 X h+ Tr[X] = S^X] . 



4 

o,/3=0 



It follows that 



l = (0|Tr[|^)(^|]|0) = i £ |(0K 1 ■ 12 

(o,/3)6Li6 



Consider a given special quadruple Q G L 16 , then we can split (6.4.1) into 
two terms: 

1 = (0|Tr[|V)(VlM 



4 E <*i**i- 



(a,/3)eLi 6 

J E k^i^i'+i E imv-'-" 2 



4 ^ 4 

(o,/3)6Q (a,/3)eL 16 \Q 



Therefore, if |?/>) and \4>) satisfy (6.34) for a given special quadruple Q, then 
(01 Co/3 IVO = f° r a ^ ( a > £0 n °t belonging to Q, i.e. the last term in the above 
equality vanishes, on the other hand K^la^l^)! 2 = 1 for all (a, (3) G Q. 
Consider now a set of positive coefficients \ a p satisfying 

E V IOK/^>| 2 < i v|^),|0)gc 4 . 

If Q C / is a special quadruple and | ■?/>), |0) satisfy ( |6.34 ), then 

i> E \(<f)\o- a fs\ip}\ 2 = E • 

(a,/3)e/ (a,/3)eQ 

□ 
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Proof: Given any four positive coefficients {A Q/ g} (a,/3)ei satisfying (6.24) 
from the assumption and the previous Lemma it follows that: 

(a,/3)Gl=Q aQ r, 

since J is a special quadruple for each of its points (a ,f3 ). Summing 
over all (a , 0o) 6 !■> eac h point is counted four times, thus 

4 = Yl J] ^<4 

(a,/3)£/ (ao,^o)G/ (a,/3)eQ ao/3o 



whence the criterion (6.24) for separability is satisfied. □ 



6.4.2 Separability and Quadruples 

The main tool in the previous proof is the fact that all points in I belongs 
to a same number, 1, of special quadruples, / itself; that is each point in / 
belongs to a special quadruple contained in /. 

Because of its importance in what follows, we now study this occurrence 
in more detail considering a higher rank lattice state corresponding to the 
subset 

/ = {(0, 0), (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3)} 
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3 



(6.41) 
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For (0,0), of the set (6.40) the following quadruples are contained in /: 
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X 


1 
1 




X 
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3 



(6.42) 



Given another point / 3 (a, (3) ^ (0,0), consider the lattice state pi„ = 
a a/ 3 (g) tpi<T a/ 3 (g) 1, where J Q/ g C L 16 is the set / transformed into r a p[I] by the 
map r a p introduced in Lemma ( |6.40 ). The point (a, (3) G / is mapped into 
(0, 0) G Jq./3 and, by Lemma Q, it belongs to 6 special quadruples contained 
in I a/ 3 that are images under r aj 3 of the 6 special quadruples relative to 
(0, 0) G /. Therefore, each point (a, (3) E I belongs to 6 special quadruples 
contained in /. 



Example 36 As a concrete example, consider the following quadruples 
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(6.43) 

They are contained in I but do not belong to Q 00 ; consider the point (3, 3) G 
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Qi, it is mapped into (0,0) by r 33 which sends I and Qi into 



h3 = T~3 3 [I] = 
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733 [Qi] = 
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One thus checks that t 33 [Qi] is a special quadruple contained in the trans- 
formed subset I 33 . Analogously, t 32 , respectively r 12 , map (3,2), respectively 
(1,2), into (0,0) and I , Q 2 , respectively Q%, into 



I32 = r S2 [I] = 



T32IQ2} = 



hi = r 12 [I] = 



ri2 [Q 3 ] = 
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G Q, 
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e Q, 
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Definition 29 Given a subset I C L 16 we shall term a covering of I any 
collection Qi of special quadruples (not necessarily disjoint) contained in I 
and denote by Nq t its cardinality. Further, we shall denote by the 
number of special quadruples in Qi that contains the point (a, (3) £ I and 
refer to them as to the multiplicities of Qj. Finally, we shall call uniform any 
covering Qi of I of constant multiplicity: M^g = Mq t , for all (a,/3) G /. 



The usefulness of uniform coverings can be seen as follows: summations 
over (a, (5) G / can be split into sums of contributions from the special 
quadruples of any covering by taking into account to how many special 
quadruples a point (a, (3) does belong: 

Y M S ^2 Y ' 

{a,0)el QeQi(a,/3)eQ 

whence, if the covering is uniform with multiplicity Mq I7 

A ^ = mT a ^ • 

{a,0)&I 21 QeQi (a,/3)gQ 



Lemma 7 Let I C L\ G contain Ni points and Qj be a uniform covering of 
I of cardinality Nqj and multiplicity M Ql ; then, 



N c 



M Ql Nj 
4 



(6.44) 



Proof: Each of the Nj points in I belongs to ni special quadruples each 
containing 4 points of I. □ 



In the case of (6.41), the special quadruples in (6.42) form a covering of 



/, but not a uniform one as its multiplicities are 



n 



oo 

Qi 

ii 



71 



33 
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n 



Qi 

22 

Qi 
21 



n. 



n 



Qi 
23 

Qi 

12 



n 



32 



n 



Qi 
13 



n 



Qi 
31 



A uniform covering I is provided by 2 special quadruples in (6.42) and the 
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3 in (6.43): 
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(6.45) 



Since all points of / belongs to exactly 2 special quadruples contained in /, 
this uniform covering has multiplicity nf = 2: 

4 x N Ql = 20 = mNj = 2 x 10 . 



This is also a minimal covering among the uniform ones; indeed, by (6.44) 
as Nq t is an integer, 4Nq = 10 x m can only be satisfied by rij 
N 3 m > 1. 



2m, 



It thus follows that the lattice state corresponding to the subset (6.41) 
is separable as it can be decomposed into a convex combination of PPT 
separable rank-4 lattice states: 

1 1 5 

Pi = ^(^00 + ^11 + ^12 + ^13 + ^21 + ^22 + ^23 + ^31 + ^32 + ^33) = 5 ' 



where {Qj}^ =l are the special quadruples of the minimal covering (6.45). 

The previous argument proves the separability of the lattice state (6.41), 
an issue which left unsettled in [5]; it can be generalized as follows. 

Proposition 6 Suppose I C L\q is a subset with minimal covering Qi of 
cardinality Nq { and multiplicity Mq i; then pi can be convexly decomposed as 

N Qj 



Pi 



1 



3=1 



where Qj C I are the special quadruples in Q, and is thus separable. 
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Proof: Let Qj, 1 < i < Nq i} be the elements of the minimal covering Qj. 
From Lemma (|6| it follows that 

N Ql 

^ ^ E E x ^ = m qi E 

3=1 (a,/9)6Qj (a,p)el 

N Ql _ Ni 
M Ql 4 



E A 

(o,/9)e/ 



a/3 



< 



so that Corollary ^ ensures separability. Furthermore, using (6.44), 

pi = Wi E p ^ 



4 



^7 j 

M Q/ iV Ql E 4 E P ^ 

M 21 3=1 (a,/3)eQj 



1 S/ 

r~ E ^ • 

i j= i 



□ 



Example 37 The following lattice states mentioned in f^j were also un- 
known to be separable: 
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;(P00 + Pll + Pl2 + Pl3 + P23 + P31 + P32 + P 
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The next state is: 

P9 = l(Poo + P11 + P12 + P13 + P21 + P23 + P31 + P32 + P33), 
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From Proposition ^ one draws the following conclusions: 
Corollary 4 All the lattice states with Nj > 14 are separable. 
Proof: 

1. iV/=16: is trivial. 

2. iV/=15: As we have seen in (section III of this chapter), every point 
(a, (3) can take part in 15 quadruples. Therefore in Li 6 we have: 



Nr 



16 x 15 



By dropping out one point from the lattice, we decrease the number 
of quadruples by three for each point and therefore each point has 12 
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quadruples contained in I: 



Nr. 



15 x 12 
4 ' 



3. Nj = 14: Here there are two points missing from Lie, note that every 
two given points can share three quadruples. This is easy to verify it; 
by sending one of the two into the origin (0, 0), one sees that the other 
one appears in exactly three Q's (section III). Therefore, there are six 
pairs (a, (3) £ J, where the two points (a%, fii) and (a 2 , A2), remove only 
five quadruples, leaving 10 special quadruples contained in I. While for 
the other eight pairs (a, 0) £ J, they remove six quadruples. Therefore 
any pi4 would look like: 
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Then one can divide it into two subsets: I\ : the six elements with 10 
special quadruples contained in I, and J 2 : the eight elements with 9 
special quadruples contained in I: 
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In order to show that all p%& are separable, we need to show that they 



satisfy (6.23): 



< 



14 

T 



a,pei 

for all choices of the coefficients \ a p introduced in Corollarj|2j But 
since we can divide pi4 into two subsets, we have: 



a,PeI ct,peh a,p&I 2 



< 



< 



4 4 
N h + N h 



14 

T' 
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Where we have used the separability of the two states pi x and pi 2 . 
Indeed in case of Nj = 14 the minimum number of special quadruples 
for every point to be included in I, is equal to n* = 2, which is much 
less than what happens. 

□ 



6.5 Entanglement and lattice geometry 

In this section we single out a family I spec of geometric patterns such that 
if / G I S pec, then pi is surely entangled. Given a sub-set / C L 16 , choose A Q / 
such that for a certain (a , (3 ) £ I 

1 + 5 



A 



K/3 = j V(a,/3) ^ (a ,/3o) 



(6.47) 



where 5 > is a suitable parameter. Then, A Q/ 3 = — ^— — > and 



(a/3)e/ 



inequality (6.25) in Corollary |2j is satisfied. According to the same corollary, 
given such a subset J, the corresponding lattice state pi is entangled, if also 
inequality (6.24) is fulfilled, namely if 



35>0suchthat^|(^| ( T Qo/3o |V)| 2 + ^ Yl 



4' 



a. 



• a ^)| 2 <l, W\ v ),\^)eC 4 



(6.48) 



Remark 6 Suppose that inequality (6.48) can be satisfied by 5 > uniformly 
in the vectors \ip) and \<p). Then, using (4.25) and (4.26), the choice of 



coefficients in (6.47) corresponds to a witness of the form 



A = Tr — A cp = - SaP - 



8 



«o/3o 



(6.49) 



By comparison with the expression of the transposition in (4.23), one sees 



that, unlike the 6 negative contribution of the latter, A presents only one 
negative contribution. Yet, for the special subset T spec , it proves to be more 
flexible as an entanglement witness. 
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Clearly, if (a , (3 ) belongs to a special quadruple contained in /, such a 
5 > cannot exist; indeed, if \ip) and \ip) satisfy (6.34), t hen | ( y^cy,/^ [■;/>) | > 

and - | (ip\a a p\ip) | 2 = 1, together with inequality (6.48) yield 5 = 0. 

(«,/3)e/ 

We now show that lattice state pi for which not all points of / belong 
to special quadruples contained in J, such a 5 > can indeed be found and 
thus that they are entangled. 



Theorem 20 Given a lattice state pj, if there exists a point (ao, (3q) G / 
such that Q ao/ 9 ^ / for all Q Qo /3 G Q, then pi is entangled. We shall call 
any such I a special sub-set and denote by X spec their family. 



It is very constructive to present some examples before going through the 
proof. 

Example 38 According to the last theorem the following state is entangled, 
as non of the special quadruples of (0, 0) are included in I: 
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(6.50) 



Another example of such entangled lattice states is the following: 
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(6.51) 



As we observe the point (3,3) has no special quadruple contained in I, it can 
be easily verified when this point is sent to (0,0), the resulting state is: 



iVj = 10 
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As the last example consider the following state: 



Nj = 11 
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(6.52) 



Also in this case (0, 0) has no special quadruple contained in I. 



Proof of Theorem | 

According to Lemma |4j the point (o>o,/3q) can be transformed into (0,0) 
and / into a new set, that we shall denote again by I for sake of simplicity, 
without altering the entanglement or separability of the transformed pi with 
respect to the initial one. Then, the assumption of the theorem translates 
into the fact that no special quadruple in the list (6.40) is contained in /. 
Having set (ao, A)) — (00), inequality (6.48) now reads 



where Af v = - ^ \(<f\cr^ 



(6.53) 



It proves convenient to introduce the following -^-dependent 4x4 matrices 

(6.54) 



so that one has to prove that 



35 > such that A+f = {(p\A* s \(p) < 1 V 



Further, from (4.25) 



Af + A% 



< AJ < 1 



(6.55) 



(6.56) 



Clearly, the major obstruction to Aff < 1 with 5 > arises when Aj 
has eigenvalue 1. Then, because of (6.56), the corresponding eigenvectors, 
At\ ( p) = \<p)> satisfy 



cr, 



a/3 







\<P) 



(a,/3)e/ c 



V(«,/3)e/ c . 

(6.57) 
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Let the eigenvalues of A]f be decreasingly ordered and consider the spectral 
decomposition 



(6.58) 



where Pf(l) projects onto the eigenspace relative to the eigenvalue 1 and 
Pf(j) are the other orthogonal spectral projectio ns rela tive to the eigenvalues 



< dj(j) < 1. Then, if PJ (1) ^ 0, 5 > in (6.55) is only possible with 



Pj(l)\ip) = 0. These preliminary considerations al 
through a series of lemmas and corollaries. 



ow us to prove ( 6.48| 
□ 



Lemma 8 With the notation of (6.58), if Pf(l)\tp) = for all \if>) G C 



then (6.55) is satisfied. 



Proof: If Mi 



sup 



\RJ 



1, then, by compactness, there exists a 



converging sequence ip n — > ip* of vectors in C 4 such that Rf n _L P/ n (l) and 
||i?j n || — > 1. Then, A^ n converges in norm to Af with \\Rf \\ = 1 and 
Rf JL Pf{l) which is a contradiction. Therefore, Mj < 1; hence, choosing 
< 5 < 4(1 - Mj), from Pf(l)\ip) = and Pf(l)Rf = 0, one gets 



I At 



max 



S 



< max 



< 1 . 



□ 



Corollary 5 Given a lattice state pi and G C 4 ; let Vic be the subspace 
spanned by the vectors c a ^\ip) with (a, (3) G I c , where I c is the complement 
of I. IfVfc = C 4 for all \ip) G C 4 , then pi is entangled. 



Proof: From ( |6.57| ) it follows that Pf(l] 
Lemma |8] applies. 



for all 



G C 4 ; thus 
□ 



Based on the previous two results, we now focus upon when Pf(l) ^ 
and show that it projects onto a subspace orthogonal to 



Le mma 9 If a^ip) = ±\if)) for some {p,v) G Li 6 , then, with the notation 

ofEsik, pf(i)iv>) = o. 
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Proof: If (n,v) G I c and (ip\Af\ip) = 1, then, from (6.57), \(p) JL (J a p\ip) 
for all (a, /3) G J c , hence to ±|V") — ^^A^)- Suppose then that (/i, v) G / 



and rewrite inequality ( |6.55[ ) as 

■mm 2 + \ E i^i^i^)i 2 + z E k^k^i 



" ' ' ' - • > 1 ^ K ■ » 1 v * 2 < 1 

4 1 - r 1 "■'>-/> - 1 ' 



where the index set I has been split into 

h = {(a, (3) E I : [a af3 , a^) = 0} and 7 2 = {(a, 0) G I : {a a/3 , a^} = 0} . 

The vectors <J a p\i>) from these two subsets are orthogonal; indeed, [a-ys, c^v] = 
and {a a p, a^} = yield 

Therefore, the following two are orthogonal matrices: 

/i3(a^)^(00),( A t, l /) 



\ E a <*pW)^\ a <*P ■ (6.60) 



4 

(a,/3)e/ 2 



Then, A^|0) = |0) can only be due to ' \(p) = \ip), in which case \<p) J_ |^/>). 
Indeed, || Aj || < 1. 

This can be seen as follows: at most three a a p may contribute to the sum 
in Aj and they must anti-commute. In fact, if there were two commuting 
<T a /3 contributing to the sum, then, as they would commute with and 
form a special quadruple Qoo contained in I which is excluded by hypothesis. 
Therefore, the <J a p contributing to the sum must anti-commute and, from 
Lemma [5j cannot be no more than three. Suppose this is the case; denote 
by S a , a = 1,2,3, these three a a p such that {a a p, a^ u } = and rewrite 

1 1 3 

Al = -\ij)^\ + -J2s a \^)(^\S a , S = U. 

Without restriction, we choose \ip) such that cr^ip) = \ip). Each S a \ip) is an 
eigenstate of a^ u belonging to the same twice degenerate eigenvalue 1. Let 
P project onto the corresponding eigenspace; then, [S a , P] = and the rank 
2 matrices T a = P S a P = S a P = P S a satisfy the Pauli algebra (2.10) over 
C 4 . Thus, (2.11) holds on C 4 with T a replacing a a , yielding 

3 

K = + \ E s a mms a = \wm + \p . 

a=0 
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If there are less than three anti-commuting contributions S a , a^O, then the 
second equality becomes a strict inequality. Therefore, ||A| ' || < 3/4 < 1. □ 

Let us now consider the set of a a p indexed by the points (a, (5) in the 
complement set I c , that we list as (aij, 1 < i < N — Nj. The following 
facts are a consequence of the structure of Q 00 studied in Lemma [5] and 
concern the points of I c and the special quadruples in Q 00 they belong to. 

1. Take the first point (ax,Pi) £ I] then, it belongs to three quadruples 



in (6.40) which cannot be contained in I. 



2. Also the second point («2, P2) £ I c may eliminate three quadruples from 
those contained in J; it does so, only if it does not share any special 
quadruple with those of (ai,/3i): this can happen only if o ax ^ and 
a a2 f3 2 anti-commute. If they commute, then and {0.2,^2) share 

one special quadruple. In this case, (a^j P2) eliminates only two special 
quadruples from those contained in /. It thus follows that the minimum 



cardinality of I c complying with the hypothesis of Theorem 20 , namely 
that no special quadruple of (0,0) is contained in J, is five with the 
corresponding a a p forming an anti-commuting set. 

3. Consider the third point (a 3 , (3 3 ) G I c : the number of special quadruples 
eliminated by this point from those contained in / is equal to three 
minus the number of previous points whose a a p commute with cr a3 p 3 - 

4. Moving to the fourth point (04, ^4) and further on, the rule is the 
same; only, one has to take into account that u ai p i does not commute 
with more than three a Qj g. Therefore, if this happens, (a 4 /3 4 ) does not 
eliminate from those contained in / any further special quadruple with 
respect to those already eliminated by the previous points. 

Based on these properties, we now distinguish the following cases. 

Case 1. card(/ c ) = 5: the matrices {c a< ft}f=i> with (c^,/^) G I c anti- 
commute. 

Case 2. card(/ c ) > 5 and there is a set of anti-commuting {caift}f=ij with 
(ati, (3i) G I c : let (a^, f3§) G I c not belong to the index set of the five 
anti-commuting cr's . The corresponding must commute with at 
least one of them and may commute with at most three of the a ai ^. 
Therefore, there must exist at least two of them that anti-commute with 
<7 a6i 8 6 . By putting together cr a6( g 6 with the one in the anti-commuting 
set which commutes with it and with the two which anti-commute with 
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it, then one can construct a set K of four a a p where 3 of them anti- 
commute among themselves and the fourth one commutes with only 
one of them. 

Case 3. card(/ c ) > 5, but there is no set of 5 anti-commuting a ai p t indexed 
by (aiijPi) G I c . If there are 4 anti-commuting a a p indexed by points 
in I c , these latter correspond to twelve special quadruples in Q o n °t 
contained in /. Therefore, there are three special quadruples of Q 00 left 
that also can not be contained in /. Thus, some of their points must be 
excluded from /. These three quadruples share a common point that 
cannot be excluded from J, otherwise the corresponding a a p, together 
with the four which have already been assumed to anti-commute, would 
constitute a set of five anti-commuting cr's contrary to the hypothesis of 
the present case. Therefore, in order to eliminate the three remaining 
special quadruples from being contained in I, one point for each of the 
three special quadruples must be excluded from I. The corresponding 
<j a/ 3 will then anti-commute. Then, we have constructed two sets of 
it's, the first one containing four anti-commuting a a p , the second one 
three anti-commuting a a p. We already know that each element from Vi 
commutes at least with one from U\\ further, it cannot commute with 
more than two of them. Indeed, if this happened, one would have two 
anti-commuting a a p in the same special quadruple which is impossible. 

This makes it possible to contruct a set K as in the previous case con- 
sisting of three anti-commuting a a p and a fourth one which commutes 
with only one of them. This set K can be constructed as follows: con- 
sider the set U\ with a ai p t , i = 1,2,3,4, take the first three of them 
and a a ipi from U-i- if <J a 'i3' commutes with only one cr cti p i , i — 1, 2, 3, the 
set K consists of these four cr's. Otherwise, if a a >pi commutes with two 
(T ai p v then replace one of them with <J ai p 4 ; as said before a a /pi does not 
commute with it because of the previous argument. Therefore, these 
four elements form the set K sought after. 

The above procedure works also in the case of sets of anti-commuting 
cr's with less than four elements. 

The next two Lemmas concern Case 1 above, where I c corresponds to a 
set of 5 anti-commuting a a p\ the first lemma specifies the structure of the 
anti-commuting set, while the second one is about the sub-space generated 
by these matrices acting on a vector state \ip). 
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Lemma 10 All sets of 5 anti- commuting o~ ae i3 e , Pe) G I c , must be of the 
form 

Ki = {(T i , <?0j , CTlk , 0~ 2k , <?3k} Or K 2 — {<TiO , <JjO , CTfci , <7fc 2 , <T k3 } 

wift k such that the anti- symmetric tensor e^ k 7^ 0. 

Proof: Since we have 5 indices ae and (3e to choose among 0, 1, 2, 3, at least 
two a s and two (5 s must be equal: let us consider a± = a 2 . 

If a± = a 2 = 0, in order to have {ctq^ , o"o/3 2 } — 0, the corresponding 
fii and 02 must be different and different from 0; therefore, f3\ = i, (3 2 = j, 
with i 7^ j, i 7^ 0, j 7^ 0, so that the first two anti-commuting matrices are 
a 0i and a j. In the three (ae/3e) left, there cannot appear ae = 0, otherwise 
the corresponding /3g, which cannot be and must equal either % or j, would 
lead to a o at p t violating anti-commutativity. Then, with the three remaining 
an 7^ 0, the three corresponding j3g must be different from 0, % and j, otherwise 
c ae i3 e would commute with either <7oi or a j. Thus, /?3 = /?4 = /?5 = (3 k with 
k 7^ i, k 7^ j , while the corresponding ae must be different, namely o 3 = 1, 
0:4 = 2 and o 5 = 3. This gives the set K\ 

If a 1 = a 2 = k 7^ 0, the corresponding /3 12 cannot be equal and must be 
both different from 0. Thus, two (J ai p t in the anti-commuting set are a ki and 
(jfcj, with i j = 1,2,3. Only one of the remaining ae can be 0, otherwise 
we would be back to the previous case: for aop e to anti-commute with a k % 
and a k j, Pe must be equal to s with s ^ i, s ^ j . In order to anti-commute 
with a ki , a k j and cr 0s , o~ a4/ 3 4 and <7 a5j g 5 , with 0:45 7^ 0, must be of the form 
<j pi and (T ? j with p 7^ g and both different from and k. But then, they 
would commute among themselves; thus, 03,4,5 7^ 0. At least one of them, 
say o 3 must equal k; then, anti-commutativity only holds if ^3 = r 7^ 
with % 7^ r, j 7^ r. Finally, since at least two /3's must be equal, in order to 
anti-commute among themselves and with a ki) a k j and a kr , a a4 p 4 and cr a6 p 5 
must have /3 4 = /3 5 = and o 4 = m, o 5 = n such that the antisymmetric 
tensor e kmn 7^ 0. This fixes the set K 2 . The same result would obtain if one 
started arguing with two equal (3 s instead of a s. □ 

Lemma 11 Let K 12 be the two sets of 5 anti-commuting a a/ 3 as described 
in the previous lemma. For any \ip) G C 4 ; let 

Kf j2 = Linear Span{a a p\il>) : a a p G K h2 } . 
Then, any vector orthogonal to Kf 2 is also orthogonal to 
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Proof: Consider the set {a 0i , a j , a\ k , a 2k , o"3fe} and let |0) 2 and |1) 2 be 
the eigenvectors of ov Furthermore, fix % and j so that 



cr k \0) 2 


= |0> 2 , 


^11)2 = 


|1) 2 ; 


CTi\0) 2 


= |1>2, 


O r i| 1 )2 = 


|0) 2 ; 


^|0>2 


= i|l) 2 , 


^•|1)2 = 


-i|0) 2 



Given l^) , \<p) e C 4 , consider their expansions with respect to the orthonor- 
mal basis {\ij) = <g> \j)2}}j=o where |z)i, % = 0, 1, are eigenvectors of 03: 

Wl = 5Zjj=o^ij|y) and Iv 9 ) = Z)ij=o Actin g with first four a a/3 

, we have: 



on 



C0i\ 

cr j\ 

Clfcl 
0"2fc| 



^oo|01)+Voi|00)+^ 10 |ll)+^ii|10) , 
z(^oo|01> - ^oi|00) + Vio|ll) - Vn|10» 
^oo|10)-Voi|H) + ^o|00)-Vii|01) , 
i(Voo|10) - ^oi|H) - V'ioIOO) + V11IOI}) 



(6.61) 



If \tp) is orthogonal to each of the previous vectors, it follows that 

<Poi -000 + ^jT0io = , ^ooV'oi + ^To^n = 
^loV'oo - ^iT0oi = , ^00^10 - ^oiV'ii 
These relations recast in matrix equation read 



. 








"000 




-^11 



ip 10 

^11 

^00 —001 











<Poi 




<Pio 







. 



The determinant of the matrix is 2^ O oV'oiV ; io'0iij thus, \ip) orthogonal to the 
linear span of the vectors in (6.61 ) can only exist if all components of \ip) are 
non-zero. Then, 



\<P> 



10 



V?oo 



11 



<Pio 



01 



000 ' ' 000 r '" '•'on 

^|00)+^|01)-^|10)-^|1) . 



and 



Such a vector results orthogonal to both the fifth vector \ip) and 
Similar considerations hold for the set {a i0 , <Tj , a k i , o k2 , cr fc3 }. 



itself. 

□ 



The next Lemma concerns the dimensionality of the linear spans of vectors 
of the form a^ip) where the matrices belong to the sets discussed in Cases 
2 and 3 in the proof of Theorem 20 where three anti-commute and the fourth 
one commutes with only one of them. 
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Lemma 12 Consider a sub-set K = {<7 ai ^ : i = 1, ...,4} consisting of three 
anti- commuting matrices a a p plus a fourth one which commutes with only 
one of these three, say the first: 

Wa^ , (Tojb} = , i,j = 1,3 [a Q4/?4 , a ai/3l ] = . (6.62) 

Then, unless \ip) G C 4 is an eigenstate of some a^ u , Vk = C 4 , where Vk is 
the linear span of {cr Q , i/3i }f =1 . 



Proof: Suppose \tp) is not an eigenstate of any a^, then the vectors 
CaijSilV'}) a aili-M ) ) cannot be proportional. For the same reason, the vectors 
CaiySilV')) GuiP-M?) an d cr o4 J 8 4 |'0) are linearly independent. Indeed, suppose 

°<upM>) = ctcr ai iiM) + /3cr ai/9l |^) . 

Then, acting on both sides with a Q4/ 3 4 we obtain \tp) = a Si\ip) + (3S2\ip), 
where Si = o a ^^a a ^ % , % = 1,2. When substituting this expression for \ip) in 
the right hand side of the equality, the relations (6.62) yield 

|^) = ( a 2 - f3 2 )\ip) + a/3 [a a%Pl , a a2 p 2 ]\^) , 

which, for non-trivial a/3, is only possible if is an eigenstate of the the 
a^v proportional to the commutator in the previous equality. The same 
conclusions can be drawn about the linear independence of the four vectors 
spanning Vk- Indeed, if 

the same argument of before obtains 

= -( a 2 + (3 2 + -f 2 )\ip) + a^{a aiPl , a Q4/34 }|^) • 

This implies that, unless is eigenstate of the a^ v arising from the anti- 
commutator in the above expression, the four vectors spanning Vk are lin- 
early independent. □ 



The previous results show that, as a consequence of the assumption that 
no special quadruple Q G <2oo is contained in /, given a vector G C 4 
the sub-space linearly generated by (T a p\i/)) with (a/3) G I c is either C 4 or 
contains \ip). The proof of Theorem 20 is thus completed by using Lemma [8] 
and its corollary. □ 
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It might look like that, when all the points of the lattice have at least one 
special quadruple contained in I, the corresponding lattice state might be 
separable. But in the following examples, we will show that this is not in 
general true. 

Example 39 Consider the following state: 



N r 
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X 
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1 
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X 


X 







1 


2 


3 



This looks very similar to the state in example \6.5S\ where none of the special 
quadruples of (0, 0) were included in I. However in the above lattice state, 
(0, 0) has two special quadruples contained in I, which are: 
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1 
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3 
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But from Proposition^ we know that this state is NPT: the contribution to 
I from the row and the column for the point (3, 0) is 6 > y . 
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Summary and Outlook 



In this thesis work, we have studied the role of positive and completely pos- 
itive maps in detecting entanglement. In the first chapters we have reviewed 
the necessary techniques that have been employed in the last chapter. 

Indeed, in the final chapter of the thesis we have considered a particular 
class of bipartite states pi on C 4 ®C 4 , introduced in jl], called Lattice States. 
These are uniform mixtures of projections indexed by points (a, (3) belonging 
to subsets I of the finite square lattice of cardinality 16. The projections are 
generated by the action of I4 <g> <j a p on the completely symmetric vector in 
C 4 ® C 4 , where cr a p is the tensor product of the Pauli matrices <r a , ap . We 
have generalized them to states on C 2 ™ <8> C 2 ", referred to as a— diagonal 
states. 

One of the main issues was to decide whether a given a— diagonal state, 
is entangled or separable. We have tackled this problem using the results 
of [65] and showed that, starting from a general non-diagonal positive map, 
possible entanglement of a— diagonal states would be revealed using only a 
particularly simple sub-class of positive maps adapted to the diagonal struc- 
ture of the states, that are combinations of elementary maps of the form 
X 1 y a a/ 3 X a a p. 

Back to the Lattice States in C 4 <g> C 4 , we have shown the compatibility 
of the results given in [5], E] , with ours. 

As the next step, we have partially shown how the entanglement (respec- 
tively, separability) of the Lattice States can be explained through their 
geometrical structure. The results obtained are based on particular separable 
lattice states consisting of four points, called special quadruples, which play 
a crucial role in this game. 

Using this notion, we have shown that if there exists a point in the set 
/, such that none of its special quadruples are contained in J, then the 
corresponding Lattice State pi, is entangled. We have also provided a proper 
entanglement witness for such a state based on the results obtained in the 
first part of the chapter. 

Also, using the notion of uniform covering by special quadruples con- 
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tained in the subset /, we could show that some of the Lattice States whose 
separability was left unknown in previous works, are indeed separable: 
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However, there are still some Lattice States whose separability could not 
be decided by the methods developed in this thesis; for instance, 
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For this Lattice State, all points have at least three special quadruples con- 
tained in I, but no uniform covering could be found, nor could be excluded. 



Another interesting issue, still to be understood, is the relation between 
the special quadruples and partial transposition. The difficulties can be ap- 
preciated by looking at the following state: 
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which is NPT entangled, nevertheless all points belong to at least one special 
quadruples contained in /. Obviously, in this case there does not exist a 
uniform covering of special quadruples as the state is entangled. 



Further interesting issues to be tackled are whether separable Lattice 
States carry quantum discord and whether the structure of PPT entangled 
Lattice States could be investigated in terms of Edge States and Unextendible 
Product Bases. 



Appendix A 
Convexity 



A good introduction of Convex Analysis is given in [TJ [58] . 

Definition 30 Convex Set 

Let S be a finite set of points in IR d . S is said to be convex if 

ax + (1 - a)y G S Vx,y G S < a < 1. (A.l) 

Consequently a point y G M. d is said to be a convex combination of {xi}f =1 
C R d if: 

y = ^PiX h ^pi = l, ft > fori = l,...,n. (A.2) 

j i 

Example 40 The set ofstates and the subset of separable states p s G M n (C) 
are convex sets. 

Definition 31 Convex hull 

Convex hull of a convex set S, denoted by conv(S) is the set of all the 
convex combinations of the elements of S: 

conv(S) = {^piXi, y^^Pi = 1, Pi>0xieS}. 

i i 

Convex hull of S is the smallest convex set containing S. 

Let S be a convex set. Then its elements split into two subsets: 

Definition 32 Interior points 

The point x G S is an interior point if the line segment from any point 
in S to x can be extended within S . 
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Definition 33 Extreme points 

A point y G S is said to be extreme if: 

y = axi + (1 — a)x 2 , x±,x 2 G S ==>• (y = X\ A a = 1) V (y = x 2 A a = 0). 

In other words, extreme points of S do not belong to any proper segment 
defined by two distinct elements of S. 

Example 41 The pure states or projectors are extreme points for the set of 
density matrices, whereas the density states of mixed states are the interior 
points. 

When an inner product (.|.) is defined on C, we have: 

Definition 34 Hyperplane 

Let a and xq be vectors in lR d . Then a hyperplane is a set of the form 
{x\(a\(x — x )) = 0}. The vector a is called the normal vector for the hyper- 
plane. 

The hyperplane divides M. d into two halfspaces: 

Definition 35 Half space A halfspace is a set of the form: 

{x\(a\(x - x )) < 0} 

The concept of hyperplanes and halfspaces are very useful as using them 
we can separate the convex sets which do not intersects: 

Theorem 21 Let C and D be two convex sets such that C fl D = 0. Then 
there exists a / and x such that (a\(x — x )) < b for all x G C and 
{a\(x — x )) > b for all x G D. In another words, (a\(x — x )) is non-positive 
on C and non-negative on D. 
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If a hyperplane is tangent to the set, at some extreme point x , then it is 
called supporting hyperplane. Note also that the separating hyperplanes as 
well as supporting hyperplanes are not unique. 

Definition 36 Cones 

A set C C R d is a cone if 

x G C, a > ==>- ax G C. 
C is said to be a convex cone if: 

x,yeC, a, (3 > =>- ax + (3y e C. 
To any convex cone one can associate a dual cone. 

Definition 37 Dual Cone 

The dual cone C° is defined to be: 

C° = {y.(y\x)>0, VxeC}. 
The dual cone C° is always closed, even if C is open. 
Example 42 

When the cone and dual cone coincide then the cone is said to be self-dual, 
i.e. : 

(y\x) > ViGC y > 0. 



110 



APPENDIX A. CONVEXITY 



Example 43 The cone of positive semidefinite n x n matrices is self-dual. 
This also holds for the cone of completely positive maps: CV = CV° . 

Example 44 All these concepts can be extended to the Hermitian operators 
and linear maps as well. Indeed one can easily verify that the set of positive 
operators is a convex set, and that positive (resp. complete positive) maps 
are unaffected by multiplication by positive scalars. 

Let C denote the cone of (positive, completely positive) linear maps, the 
dual cone C° in this case is defined via the Jamiolkowski Isomorphism. Let 
{Eij} be an orthonormal set of unit matrices, then the Choi matrix C$, for 
the map : M n (C) — > M n (C) is: 

C (j> = ^E ij ®(t ) {E ij ). (A.3) 



Using the Hilbert- Schmidt scalar product introduced in (2.4), the dual cone 
C° is defined to be: 

C° = {*Ij: M n (C) — ► M n (C) : Tr(C^) > 0, V0 G C}. (A.4) 

Let P(C d ) be the set of positive linear maps <fi : M n (C) — > M n (C), and 
C be a closed cone in P(C d ). Recall that the set of complete positive maps is 
denoted by CV. 
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